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NUCLEAR RECEPTOR LIGANDS AND LIGAND BINDING DOMAINS 

ACKNOWLEDGMENTS 
This invention was supported in part by grants from the National Institutes 
of Health grant number 1 R01 DK43787, and 5 R01 DK 41842. The U.S. 
Government may have rights in this invention. 

CROSS-REFERENCE TO RELATED APPLICATIONS 
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INTRODUCTION 

Technical Field 

This invention relates to computational methods for designing iigands that 
bind to nuclear receptors, crystals of nuclear receptors, synthetic Iigands of 
nuclear receptors and methods of using synthetic Iigands. 

Background 

Nuclear receptors represent a superfamily of proteins that specifically bind 
a physiologically relevant small molecule, such as hormone or vitamin. As a 
result of a molecule binding to a nuclear receptor, the nuclear receptor changes the 
ability of a cell to transcribe DNA, i.e. nuclear receptors modulate the 
transcription of DNA, although they may have transcription independent actions. 
Unlike integral membrane receptors and membrane associated receptors, the 
nuclear receptors reside in either the cytoplasm or nucleus of eukaryotic cells. 
Thus, nuclear receptors comprise a class of intracellular, soluble ligand-regulated 
transcription factors. 

Nuclear receptors include receptors for glucocorticoids (GRs), androgens 
(ARs), mineralocorticoids (MRs), progestins (PRs), estrogens (ERs), thyroid 
hormones (TRs), vitamin D (VDRs), retinoids (RARs and RXRs). The so called 
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"orphan receptors" are also part of the nuclear receptor superfamily, as they are 
structurally homologous to the classic nuclear receptors, such as steroid and 
thyroid receptors. To date, ligands have not been identified with orphan receptors 
but it is likely that small molecule ligands will be discovered in the near future for 
this class of transcription factors. Generally, nuclear receptors specifically bind 
physiologically relevant small molecules with high affinity and apparent Kd's are 
commonly in the 0.01 - 20 nM range, depending on the nuclear receptor/ligand 
pair. 

Development of synthetic ligands that specifically bind to nuclear receptors 
has been largely guided by the trial and error method of drug design despite the 
importance of nuclear receptors in a myriad of physiological processes and 
medical conditions such as hypertension, inflammation, hormone dependent 
cancers (e.g. breast and prostate cancer), modulation of reproductive organ 
modulation, hyperthyroidism, hypercholesterolemia and obesity. Previously, new 
ligands specific for nuclear receptors were discovered in the absence of 
information on the three dimensional structure of a nuclear receptor with a bound 
ligand. Before the present invention, researchers were essentially discovering 
nuclear receptor ligands by probing in the dark and without the ability to visualize 
how the amino acids of a nuclear receptor held a ligand in its grasp. 

Consequently, it would be advantageous to devise methods and 
compositions for reducing the time required to discover ligands to nuclear 
receptors, synthesize such compounds and administer such compounds to 
organisms to modulate physiological processes regulated by nuclear receptors. 

SUMMARY OF THE INVENTION 
The present invention provides for crystals of nuclear receptor ligand 
binding domains with a ligand bound to the ligand binding domain (LBD). The 
crystals of the present invention provide excellent atomic resolution of the amino 
acids that interact with nuclear receptor ligand, especially thyroid receptor ligands. 
The three dimensional model of a nuclear receptor LBD with a ligand bound 
reveals a previously unknown structure for nuclear receptors and shows that the 
ligand is bound in a water inaccessible binding cavity of the ligand binding domain 
of the nuclear receptor. 
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The present invention also provides for computational methods using three 
dimensional models of nuclear receptors that are based on crystals of nuclear 
receptor LBDs. Generally, the computational method of designing a nuclear 
receptor Hgand determines which amino acid or amino acids of a nuclear receptor 
LBD interact with a chemical moiety (at least one) of the ligand using a three 
dimensional model of a crystallized protein comprising a nuclear receptor LBD 
with a bound ligand, and selecting a chemical modification (at least one) of the 
chemical moiety to produce a second chemical moiety with a structure that either 
decreases or increases an interaction between the interacting amino acid and the 
second chemical moiety compared to the interaction between the interacting amino 
acid and the corresponding chemical moiety on the natural hormone. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a diagram illustrating computational methods for designing 
ligands that interact with nuclear receptors of the nuclear receptor superfamily. 

FIG. 2 is a schematic representation of nuclear receptor structures, 
indicating regions of homology within family members and functions of the 
various domains. 

FIG. 3 shows the aligned amino acid sequences of the ligand binding 
domains of several members of the nuclear receptor superfamily. 

FIG. 4 is a ribbon drawing of the rat TR-a LBD with secondary structure 
elements labelled. The ligand (magenta) is depicted as a space-filling model. 
Alpha helices and coil conformations are yellow, beta strands are blue. 

FIG. 5 shows two cross-sections of a space-filling model of rat TR-a 
exposing the ligand (magenta) tightly packed within the receptor. 

FIG. 6 is a schematic of the ligand binding cavity. Residues which interact 
with the ligand appear approximately at the site of interaction. Hydrogen bonds 
are shown as dashed lines between the bonding partners; distances for each bond 
are listed. Non-bonded contacts are shown as radial spokes which face toward 
interacting atoms. 

FIG. 7 is the distribution of crystallographic temperature factors in the 
refined rat TR-a LBD. The distribution is represented as a color gradation 
ranging from less than 15 (dark blue) to greater than 35 (yellow-green). 
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FIG. 8 is a ribbon drawing of the rat TR-a LBD showing the c-terminal 
activation domain to ligand. Residues which comprise the c-termina! activation 
domain (Pro393-Phe405) are depicted as a stick representation. Hydrophobic 
residues, particularly Phe401 and Phe405 (blue) face inwards toward the ligand. 
Glu403 (red) projects outward into the solvent. 

FIG. 9 is an electrostatic potential surface of the rat TR-a LBD, calculated 
using GRAPH. Negative electrostatic potential is red; positive electrostatic 
potential is blue. The c-terminal activation domain forms a largely hydrophobic 
(white). The Glu403 is presented as a singular patch of negative charge (red). 

FIG. 10 is a diagram comparing agonists and antagonists for several 
nuclear receptors. 

FIG.- 11 is the synthetic scheme for preparation of TS1, TS2, TS3, TS4 
and TS5. 

FIG. 12 is the synthetic scheme for preparation of TS6 and TS7. 
FIG. 13 is the synthetic scheme for preparation of TS8. 
FIG. 14 is the synthetic scheme for preparation of TS10. 
FIG. 15 depicts the chemical structures of several TR ligands. 
FIG. 16 is a graph illustrating competition assays in which T 3 and triac 
compete with labeled T 3 for binding to human TR-a or human TR-0. 

FIG. 17 depicts a Scatchard analysis of labelled T 3 binding to TR-a and 

TR-0. 

FIG. 18 is a chart showing the effect of TS-10 on the transcriptional 
regulation of the DR4-ALP reporter gene in the presence or absence of T3 as 
assayed in TRAFal reporter cells. 

FIG. 19 is a chart showing the effect of TS-10 on the transcriptional 
regulation of the DR4-ALP reporter gene in the presence or absence of T3 as 
assayed in TRAF/31 reporter cells. 

FIG. 20 is a chart showing the effect of TS-10 on the transcriptional 
regulation of the DR4-ALP reporter gene in the presence or absence of T3 as 
assayed in HepG2, a liver reporter cell line. 
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FIG. 21 is a partial ribbon drawing of TR-a LBD with T3 in the ligand 
binding cavity. Selected interacting amino acids are labelled, including Ile221, 
Ile222 and Ser260, Ala263, Ile299 and Leu 276. 

FIG. 22 is a partial ribbon drawing of TR-a LBD with T3 and Dimit 
superimposed in the ligand binding cavity. Interactions with Ile221, Ile222, 
Ala260, lie 299 and Leu276 are labelled. 

FIG. 23 is a partial ribbon drawing of TR-a LBD with T3, illustrating the 
three Arginine residues (Arg228, Arg262 and Arg 266 (dark stick figures)) of the 
polar pocket, three water molecules HOH502, HOH503 and HOH504, with 
hydrogen bonds indicated by dotted lines. 

FIG. 24 is a partial ribbon drawing of TR-a LBD with triac, illustrating 
the three Arginine residues (dark stick figures) of the polar pocket, water 
molecules (HOH503, HOH504 and HOH600), with hydrogen bonds indicated by 
dotted lines. 

FIG. 25 is a partial ribbon drawing of the TR-a LBD with T3 and triac 
superimposed in the ligand binding cavity. The drawing shows several interacting 
amino acid residues in the polar pocket that remain unchanged whether T3 or triac 
occupies the ligand binding cavity: Arg262, Asnl79, HOH503 and HOH504, and 
Ser277. Both Arg228 and Arg 266 occupy two different positions, depending on 
whether T3 or triac is bound. 

FIG. 26 is a stereochemical representation of the TRa LBD with Dimit 

bound. 

APPENDIX 1 is an appendix of references. 

APPENDIX 2 is a chart of amino acids that interact with a TR ligand, for 
TR complexed with Dimit, Triac, IpBr2, and T3. 

DETAILED DESCRIPTION OF THE INVENTION 

Introduction 

The present invention provides new methods, particularly computational 
methods, and compositions for the generation of nuclear receptor synthetic ligands 
based on the three dimensional structure of nuclear receptors, particularly the 
thyroid receptor (herein referred to as M TR"). Previously, the lack of three 



WO 97/21993 PCT/US96/20778 

6. 

dimensional structural information about the ligand binding domain of a nuclear 
receptor thwarted the field of nuclear receptor drug discover)', especially the 
absence of three dimensional structural information relating to a nuclear receptor 
with a ligand bound. 

Described herein for the first time are crystals and three dimensional 
structural information from a nuclear receptor's ligand binding domain (LBD) with 
a ligand bound. Such crystals offer superior resolution at the atomic level and the 
ability to visualize the coordination of nuclear receptor ligands by amino acids that 
comprise the LBD. The present invention also provides computational methods 
for designing nuclear receptor synthetic ligands using such crystal and three 
dimensional structural information to generate synthetic ligands that modulate the 
conformational changes of a nuclear receptor's LBD. Such synthetic ligands can 
be designed using the computational methods described herein and shown, in pan, 
in FIG. 1. These computational methods are particularly useful in designing an 
antagonist or partial agonist to a nuclear receptor, wherein the antagonist or partial 
agonist has an extended moiety that prevents any one of a number of ligand- 
induced molecular events that alter the receptor's influence on the regulation of 
gene expression, such as preventing the normal coordination of the activation 
domain observed for a naturally occurring ligand or other ligands that mimic the 
naturally occurring ligand, such as an agonist. As described herein, synthetic 
ligands of nuclear receptors will be useful in modulating nuclear receptor activity 
in a variety of medical conditions. 

Applicability to Nuclear receptors 

The present invention, particularly the computational methods, can be used 
to design drugs for a variety of nuclear receptors, such as receptors for 
glucocorticoids (GRs), androgens (ARs), mineralocorticoids (MRs), progestins 
(PRs), estrogens (ERs), thyroid hormones (TRs), vitamin D (VDRs), retinoid 
(RARs and RXRs) and peroxisomal proliferators (PPAP). The present invention 
can also be applied to the "orphan receptors," as they are structurally homologous 
in terms of modular domains and primary structure to classic nuclear receptors, 
such as steroid and thyroid receptors. The amino acid homologies of orphan 
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receptors with other nuclear receptors ranges from very low (< 15%) to in the 
range of 35% when compared to rat RARa and human TR-/3 receptors, for 
example. In addition, as is revealed by the X-ray crystallographic structure of the 
TR and structural analysis disclosed herein, the overall folding of liganded 
superfamily members is likely to be similar. Although ligands have not been 
identified with orphan receptors, once such ligands are identified one skilled in the 
art will be able to apply the present invention to the design and use of such 
ligands, as their overall structural modular motif will be similar to other nuclear 
receptors described herein. 

Modular Functional Domains Of Nuclear receptors 

The present invention will usually be applicable to all nuclear receptors, as 
discussed herein, in part, to the patterns of nuclear receptor activation, structure 
and modulation that have emerged as a consequence of determining the three 
dimensional structures of nuclear receptors with different ligands bound, notably 
the three dimensional structures or crystallized protein structure of the ligand 
binding domains for TR-or and TR-0. Proteins of the nuclear receptor superfamily 
display substantial regions of amino acid homology, as described herein and 
known in the art see FIG, 2. Members of this family display an overall structural 
motif of three modular domains (which is similar to the TR three modular domain 
motif): 

1) a variable amino-terminal domain; 

2) a highly conserved DNA-binding domain (DBD); and 

3) a less conserved carboxyl-terminal ligand-binding domain (LBD). 
The modularity of this superfamily permits different domains of each protein to 
separately accomplish different functions, although the domains can influence each 
other. The separate function of a domain is usually preserved when a particular 
domain is isolated from the remainder of the protein. Using conventional protein 
chemistry techniques a modular domain can sometimes be separated from the 
parent protein. Using conventional molecular biology techniques each domain can 
usually be separately expressed with its original function intact or chimerics of two 
different nuclear receptors can be constructed, wherein the chimerics retain the 
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properties of the individual functional domains of the respective nuclear receptors 
from which the chimerics were generated. 

FIG. 2 provides a schematic representation of family member structures, 
indicating regions of homology within family members and functions of the 
various domains. 

Amino Terminal Domain 

The amino terminal domain is the least conserved of the three domains and 
varies markedly in size among nuclear receptor superfamily members. For 
example, this domain contains 24 amino acids in the VDR and 603 amino acids in 
the MR. This domain is involved in transcriptional activation and in some cases 
its uniqueness may dictate selective receptor-DNA binding and activation of target 
genes by specific receptor isoforms. This domain can display synergistic and 
antagonistic interactions with the domains of the LBD. For example, studies with 
mutated and/or deleted receptors show positive cooperativity of the amino and 
carboxy terminal domains. In some cases, deletion of either of these domains will 
abolish the receptor's transcriptional activation functions. 

DNA-Binding Domain 

The DBD is the most conserved structure in the nuclear receptor 
superfamily. It usually contains about 70 amino acids that fold into two zinc 
finger motifs, wherein a zinc ion coordinates four cysteines. DBDs contain two 
perpendicularly oriented a-helixes that extend from the base of the first and second 
zinc fingers. The two zinc fingers function in concert along with non-zinc finger 
residues to direct nuclear receptors to specific target sites on DNA and to align 
receptor homodimer or heterodimer interfaces. Various amino acids in DBD 
influence spacing between two half sites (usually comprised of six nucleotides) for 
receptor dimer binding. For example, GR subfamily and ER homodimers bind to 
half-sites spaced by three nucleotides and oriented as palindromes. The optimal 
spacings facilitate cooperative interactions between DBDs, and D box residues are 
pan of the dimerization interface. Other regions of the DBD facilitate DNA- 
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protein and protein-protein interactions required for RXR homodimerization and 
heierodimerization on direct repeat elements. 

The LBD may influence the DNA binding of the DBD, and the influence 
can also be regulated by ligand binding. For example, TR ligand binding 
influences the degree to which a TR binds to DNA as a monomer or dimer. Such 
dimerization also depends on the spacing and orientation of the DNA half sites. 

The nuclear receptor superfamily has been subdivided into two subfamilies: 
1) GR (GR, AR, MR and PR) and 2) TR (TR, VDR, RAR, RXR, and most 
orphan receptors) on the basis of DBD structures, interactions with heat shock 
proteins (hsp), and ability to form heterodimers. GR subgroup members are 
tightly bound by hsp in the absence of ligand, dimerize following ligand binding 
and dissociation of hsp, and show homology in the DNA half sites to which they 
bind. These half sites also tend to be arranged as palindromes. TR subgroup 
members tend to be bound to DNA or other chromatin molecules when 
unliganded, can bind to DNA as monomers and dimers, but tend to form 
heterodimers, and bind DNA elements with a variety of orientations and spacings 
of the half sites, and also show homology with respect to the nucleotide sequences 
of the half sites. ER does not belong to either subfamily, since it resembles the 
GR subfamily in hsp interactions, and the TR subfamily in nuclear localization and 
DNA-binding properties. 

Ligand Binding Domain 

The LBD is the second most highly conserved domain in these receptors. 
Whereas integrity of several different LBD sub-domains is important for ligand 
binding, truncated molecules containing only the LBD retain normal ligand-binding 
activity. This domain also participates in other functions, including dimerization, 
nuclear translocation and transcriptional activation, as described herein. 
Importantly, this domain binds the ligand and undergoes ligand-induced 
conformational changes as detailed herein. 

Most members of the superfamily, including orphan receptors, possess at 
least two transcription activation subdomains, one of which is constitutive and 
resides in the amino terminal domain (AF-1), and the other of which (AF-2 (also 
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referenced as TAU 4)) resides in the iigand-binding domain whose activity is 
regulated by binding of an agonist ligand. The function of AF-2 requires an 
activation domain (also called transactivation domain) that is highly conserved 
among the receptor superfamily (approximately amino acids 1005 to 1022). Most 
LBDs contain an activation domain. Some mutations in this domain abolish AF-2 
function, but leave ligand binding and other functions unaffected. Ligand binding 
allows the activation domain to serve as an interaction site for essential co- 
activator proteins that function to stimulate (or in some cases, inhibit) 
transcription. 

The carboxy-terminal activation subdomain, as described herein is in close 
three dimensional proximity in the LBD to the ligand, so as to allow for ligands 
bound to the LBD to coordinate (or interact) with amino acid(s) in the activation 
subdomain. As described herein, the LBD of a nuclear receptor can be expressed, 
crystallized, its three dimensional structure determined with a ligand bound (either 
using crystal data from the same receptor or a different receptor or a combination 
thereof), and computational methods used to design ligands to its LBD, 
particularly ligands that contain an extension moiety that coordinates the activation 
domain of the nuclear receptor. 

Once a computationally designed ligand (CDL) is synthesized as described 
herein and known in the art, it can be tested using assays to establish its activity as 
an agonist, partial agonist or antagonist, and affinity, as described herein. After 
such testing, the CDLs can be further refined by generating LBD crystals with a 
CDL bound to the LBD. The structure of the CDL can then be further refined 
using the chemical modification methods described herein for three dimensional 
models to improve the activity or affinity of the CDL and make second generation 
CDLs with improved properties, such as that of a super agonist or antagonist 
described herein. 

Nuclear Receptor Isoforms 

The present invention also is applicable to generating new synthetic ligands 
to distinguish nuclear receptor isoforms. As described herein, CDLs can be 
generated that distinguish between isoforms, thereby allowing the generation of 
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either tissue specific or function specific synthetic ligands. For instance, GR 
subfamily members have usually one receptor encoded by a single gene, with the 
exception that there are two PR isoforms, A and B, translated from the same 
mRNA by alternate initiation from different AUG codons. This method is 
5 especially applicable to the TR subfamily which usually has several receptors that 

are encoded by two (TR) or three (RAR, RXR, and PPAR) genes or have alternate 
RNA splicing and such an example for TR is described herein. 

Nuclear Receptor Crystals 

10 The invention provides for crystals made from nuclear receptor ligand 

binding domains with the ligand bound to the receptor. As exemplified in the 
Examples, TRs are crystallized with a ligand bound to it. Crystals are made from 
purified nuclear receptor LBDs that are usually expressed by a cell culture, 
preferably E. coli. Preferably, different crystals (co-crystals) for the same nuclear 

15 receptor are separately made using different ligands, such as a naturally occurring 

ligand and at least one bromo- or iodo- substituted synthetic ligand that acts as an 
analog or antagonist of the naturally occurring ligand. Such bromo- and iodo- 
substitutions act as heavy atom substitutions in nuclear receptor ligands and 
crystals of nuclear receptor proteins. This method has the advantage for phasing 

20 of the crystal in that it bypasses the need for obtaining traditional heavy metal 

derivatives. After the three dimensional structure is determined for the nuclear 
receptor LBD with its ligand bound, the three dimensional structure can be used in 
computational methods to design a synthetic ligand for the nuclear receptor and 
further activity structure relationships can be determined through routine testing 

25 using the assays described herein and known in the an. 

Expression and Purification of other Nuclear Receptor LBD Structures 

High level expression of nuclear receptor LBDs can be obtained by the 
techniques described herein as well as others described in the literature. High 
30 level expression in E. coli of ligand binding domains of TR and other nuclear 

receptors, including members of the steroid/thyroid receptor superfamily, such as 
the estrogen (ER), androgen (AR), mineralocorticoid (MR), progesterone (PR), 
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RAR, RXR and vitamin D (VDR) receptors can also be achieved. Yeast and other 
eukaryotic expression systems can be used with nuclear receptors that bind heat 
shock proteins as these nuclear receptors are generally more difficult to express in 
bacteria, with the exception of ER, which can be expressed in bacteria. 
Representative nuclear receptors or their ligand binding domains have been cloned 
and sequenced: human RAR-a, human RAR-y, human RXR-a, human RXR-/3, 
human PPAR-a, human PPAR-/3, human PPAR- T , human VDR, human ER (as 
described in Seielstad et al., Molecular Endocrinology, vol 9:647-658 (1995), 
incorporated herein by reference), human GR, human PR, human MR, and human 
AR. The ligand binding domain of each of these nuclear receptors has been 
identified and is shown in FIG. 3. Using the information in FIG. 3 in 
conjunction with the methods described herein and known in the art, one of 
ordinary skill in the art could express and purify LBDs of any of the nuclear 
receptors, including those illustrated in FIG. 3, bind it to an appropriate ligand, 
and crystallize the nuclear receptor's LBD with a bound ligand. 

FIG. 3 is an alignment of several members of the steroid/thyroid hormone 
receptor superfamily that indicates the amino acids to be included in a suitable 
expression vector. 

Extracts of expressing cells are a suitable source of receptor for purification 
and preparation of crystals of the chosen receptor. To obtain such expression, a 
vector is constructed in a manner similar to that employed for expression of the rat 
TR alpha (Apriietti et al. Protein Expression and Purification, 6:368-370 (1995), 
herein incorporated by reference)! The nucleotides encoding the amino acids 
encompassing the ligand binding domain of the receptor to be expressed, for 
example the estrogen receptor ligand-binding domain (hER-LBD) (corresponding 
to R at position 725 to L at position 1025 as standardly aligned as shown in the 
FIG. 3), are inserted into an expression vector such as the one employed by 
Apriietti et al (1995). For the purposes of obtaining material that will yield good 
crystals it is preferable to include at least the amino acids corresponding to human 
TR-0 positions 725 to 1025. Stretches of adjacent amino acid sequences may be 
included if more structural information is desired. Thus, an expression vector for 
the human estrogen receptor can be made by inserting nucleotides encoding amino 
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acids from position 700 to the c-terminus at position 1071. Such a vector gives 
high yield of receptor in E. coli that can bind hormone (Seielstad et al. Molecular 
Endocrinology Vol 9:647-658 (1995)). However, the c-terminal region beyond 
position 1025 is subject to variable proteolysis and can advantageously be excluded 
from the construct, this technique of avoiding variable proteolysis can also be 
applied to other nuclear receptors. 

TR-a And TR-& As Examples of Nuclear receptor LED Structure and Function 
TR Expression, Purification And Crystallization 

As an example of nuclear receptor structure of the ligand binding domain 
the a- and 0- isoforms of TR are crystallized from proteins expressed from 
expression constructs, preferably constructs that can be expressed in E. coli. 
Other expression systems, such as yeast or other eukaryotic expression systems 
can be used. For the TR, the LBD can be expressed without any portion of the 
DBD or amino-terminal domain. Portions of the DBD or amino-terminus can be 
included if further structural information with amino acids adjacent the LBD is 
desired. Generally, for the TR the LBD used for crystals will be less than 300 
amino acids in length. Preferably, the TR LBD will be at least 150 amino acids in 
length, more preferably at least 200 amino acids in length, and most preferably at 
least 250 amino acids in length. For example the LBD used for crystallization can 
comprise amino acids spanning from Met 122 to Val 410 of the rat TR-a, Glu 202 
to Asp 461 of the human TR-/3. 

Typically TR LBDs are purified to homogeneity for crystallization. Purity 
of TR LBDs is measured with SDS-PAGE, mass spectrometry and hydrophobic 
HPLC. The purified TR for crystallization should be at least 97.5 % pure or 
97.5%, preferably at least 99.0% pure or 99.0% pure, more preferably at least 
99.5% pure or 99.5% pure. 

Initially purification of the unliganded receptor can be obtained by 
conventional techniques, such as hydrophobic interaction chromatography (HPLC), 
ion exchange chromatography (HPLC), and heparin affinity chromatography. 

To achieve higher purification for improved crystals of nuclear receptors, 
especially the TR subfamily and TR, it will be desirable to ligand shift purify the 
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nuclear receptor using a column that separates the receptor according to charge, 
such as an ion exchange or hydrophobic interaction column, and then bind the 
eluted receptor with a ligand, especially an agonist. The ligand induces a change 
in the receptor's surface charge such that when re-chromatographed on the same 
column, the receptor then elutes at the position of the liganded receptor are 
removed by the original column run with the unliganded receptor. Usually 
saturating concentrations of ligand are used in the column and the protein can be 
preincubated with the ligand prior to passing it over the column. The structural 
studies detailed herein indicate the general applicability of this technique for 
obtaining super-pure nuclear receptor LBDs for crystallization. 

More recently developed methods involve engineering a "tag" such as with 
histidine placed on the end of the protein, such as on the amino terminus, and then 
using a nickle chelation column for purification, Janknecht R., Proc. Natl. Acad. 
Sci. USA Vol 88:8972-8976 (1991) incorporated by reference. 

To determine the three dimensional structure of a TR LBD, or a LBD from 
another member of the nuclear receptor superfamily, it is desirable to co-crystalize 
the LBD with a corresponding LBD ligand. In the case of TR LBD, it is 
preferable to separately co-crystalize it with ligands such as T3, IpBr and Dimit 
that differ in the heavy atoms which they contain. Other TR ligands such as those 
encompassed by Formula 1 described herein and known in the prior art, can also 
be used for the generation of co-crystals of TR LBD and TR ligands. Of the 
compounds encompassed by Formula 1 it is generally desirable to use at least one 
ligand that has at least one bromo- or iodo- substitution at the R 3 , R 5 , R 3 ' or R 5 * 
position, preferably such compounds will be have at least two such substitutions 
and. more preferably at least 3 such substitutions. As described herein, such 
substitutions are advantageously used as heavy atoms to help solve the phase 
problem for the three dimensional structure of the TR LBD and can be used as a 
generalized method of phasing using a halogen (e.g. I or Br) substituted ligand, 
especially for nuclear receptors. 

Typically purified LBD, such as TR LBD, is equilibrated at a saturating 
concentration of ligand at a temperature that preserves the integrity of the protein. 
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Ligand equilibration can be established between 2 and 37° C, although the 
receptor tends to be more stable in the 2-20° C range. 

Preferably crystals are made with the hanging drop methods detailed 
herein. Regulated temperature control is desirable to improve crystal stability and 
5 quality. Temperatures between 4 and 25 °C are generally used and it is often 

preferable to test crystallization over a range of temperatures. In the case of TR it 
is preferable to use crystallization temperatures from 18 to 25°C, more preferably 
20 to 23°C, and most preferably 22°C. 

Complexes of the TR-a LBD with a variety of agonists, including 3,5,3'- 

10 triiodothyronine (T 3 ), 3Msopropyl-3,5-dibromothyronine (IpBr 2 ), 3'-isopropyl-3.5- 

dimethylthyronine (Dimit), and 3,5,3'-triiodothyroacetic acid (triac), are prepared 
with by methods described herein. Cocrystals of the rTR-a LBD, with ligand 
prcbound, are prepared by vapor diffusion at ambient temperature from 15% 2- 
methyl-2,4-pentanediol (MPD). The crystals are radiation sensitive, and require 

15 freezing to measure complete diffraction data. On a rotating anode X-ray source, 

the crystals diffract to ~3A; synchrotron radiation extends the resolution limit 
significantly, to as high as 2.0A for T 3 cocrystals. The composition of the thyroid 
hormone, combined with the ability to prepare and cocrystallize the receptor 
complexed with a variety of analogs, permitted the unusual phasing strategy. This 

20 phasing strategy can be applied to the ligands of the nuclear receptors described 

therein by generating I and Br substitutions of such ligands. In this strategy, 
cocrystals of the TR LBD containing four hormone analogs that differ at the 3,5, 
and 3' positions (T 3 , IpBr 2 , Dimit, and triac) provided isomorphous derivatives. 
For this set of analogs, the halogen substituents (2Br and 31 atoms) function as 

25 heavy atoms, while the Dimit cocrystal (3 alkyl groups) acts as the parent. The 

initial 2. 5 A multiple isomorphous replacement/anomalous scattering/density 
modified electron density map allowed the LBD to be traced from skeletons 
created in the molecular graphics program 05 (Jones, T. A. et al. ACTA Cryst, 
47:110-119 (1991), incorporated by reference herein). A model of the LBD was 

30 built in four fragments, Argl57-Glyl84, Trpl86-Glyl97, Serl99-Pro205, and 

Val210-Phe405, and refined in XPLOR using positional refinement and simulated 
annealing protocols. Missing residues were built with the aid of difference 
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density. The final model was refined to R^, = 21.8% and R frce = 24.4% for 
data fromtf.O to 2.2A, see Table 3. 

This phasing strategy can be applied to the ligands of the nuclear receptors 
described herein by generating I and Br substitutions of such ligands. 

Three Dimensional Structure of TR LBD 
Architecture of TR LBD 

As an example of the three dimensional structure of a nuclear receptor, the 
folding of the TR-a, LBD is shown in FIG. 4. The TR-a LBD consists of a single 
structural domain packed in three layers, composed of twelve a-helices, Hl-12, 
and four short 0-strands, Sl-4, forming a mixed 0-sheet. The buried hormone and 
three antiparallel a-helices, H5-6, H9, and H10, form the central layer of the 
domain, as shown in FIG. 4. Hi, H2, H3 and SI form one face of the LBD, 
with the opposite face formed by H7, H8, Hll, and H12. The first 35 amino 
acids of the N-terminus (Metl22-Glnl56) are not visible in the electron density 
maps. The three dimensional structure of the heterodimeric RXR:TR DNA- 
binding domains bound to DNA, amino acids Met 122 - Glnl51 of the TR DBD 
make extensive contacts with the minor groove of the DNA8. The five disordered 
amino acids (Argl52-Glnl56), which reside between the last visible residue of the 
TR DBD and the first visible residue of the LBD likely represent the effective 
"hinge" linking the LBD and the DBD in the intact receptor. 

The predominantly helical composition and the layered arrangement of 
secondary structure is identical to that of the unliganded hRXRa, confirming the 
existence of a common nuclear receptor fold between two nuclear receptors. 

The TR LBD is visible beginning at Argl57, and continues in an extended 
coil conformation to the start of HI. A turn of a-helix, H2, covers the hormone 
binding cavity, immediately followed by short 0-strand, SI, which forms the edge 
of the mixed /3-sheet, parallel to S4, the outermost of the three antiparallel strands. 
The chain is mostly irregular until H3 begins, antiparallel to HI. H3 bends at 
Ile221 and IIe222, residues which contact the ligand. The chain turns almost 90' 
at the end of H3 to form an incomplete a-helix, H4. The first buried core helix, 
H5-6, follows, its axis altered by a kink near the ligand at Gly 253. The helix is 
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composed of mostly hydrophobic sidechains interrupted by two striking exceptions: 
Arg262 is solvent inaccessible and interacts with the ligand carboxylate (1- 
substituent), and Glu256 meets Arg329 from H9 and Arg375 from HI 1 in a polar 
invagination. H5-6 terminates in a short 0-strand, S2, of the four strand mixed 
5 sheet. S3 and S4 are joined through a left-handed turn, and further linked by a 

salt bridge between Lys284 and Asp272. Following S4, H7 and H8 form an L, 
stabilized by a salt bridge between Lys268 and Asp277. The turn between H7 and 
H8 adopts an unusual conformation, a result of interaction with ligand and its 
glycine rich sequence. H9 is the second core helix, antiparallel to the neighboring 

10 H5-6. Again, two buried polar sidechains are found, Glu315 and Gln320. 

Glu315 forms a buried salt bridge with His358 and Arg356. The oxygen of 
Gln320 forms a hydrogen bond with the buried sidechain of His 175. The chain 
then switches back again to form H10, also antiparallel to H9. Hll extends 
diagonally across the full length of the molecule. Immediately after Hll, the chain 

15 forms a type II turn, at approximately 90* to Hll. The chain then turns again to 

form H 12, which packs loosely against H3 and HI 1 as part of the hormone or 
ligand binding cavity. The final five amino acids at the C-terminus, Glu406 - 
Val410, are disordered. 

20 TR LBD's Ligand Binding Cavity As An Example Of A Nuclear Receptor's 

Buried Ligand Cavity 

The three dimensional structure of the TR LBD leads to the startling 
finding that ligand binding cavity of the LBD is solvent inaccessible when a T3 or 
its isostere is bound to the LBD. This surprising result leads to a new model of 

25 nuclear receptor three dimensional structure and function, as further described 

herein, particularly in the sections elucidating the computational methods of ligand 
design and the application of such methods to designing nuclear receptor synthetic 
ligands that contain extended positions that prevent normal activation of the 
activation domain. 

30 Dimit, the ligand bound to the receptor, is an isostere of T 3 and a thyroid 

hormone agonist. Therefore the binding of Dimit should reflect that of T 3 , and the 
Dimit-bound receptor is expected to be the active conformation of TR. The ligand 
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is buried within the receptor, providing the hydrophobic core for a subdomain of 
the protein, as shown in FIG. 5 a and b. H5-6 and H9 comprise the hydrophobic 
core for the rest of the receptor. 

An extensive binding cavity is constructed from several structural elements. 
The cavity is enclosed from above by H5-6 (Met 256- Arg266), from below by H7 
and H8 and the intervening loop (Leu287- Ile299), and along the sides by H2 
(185-187), by the turn between S3 and S4 (Leu276-Ser277), by H3 (Phe215- 
Arg228), by Hll (His381-Met388) and by H12 (Phe401-Phe405). The volume of 
the cavity defined by these elements, calculated by GRASP (Columbia University, 
USA) (600 A3), is essentially the volume of the hormone (530 A). The remaining 
volume is occupied by water molecules surrounding the amino-propionic acid 
substituent. FIG, 6 depicts various contacts (or interactions) between TR's LBD 
and the ligand. 

The planes of the inner and outer (prime ring) rings of the ligand are 
rotated from planarity about 60° with respect to each other, adopting the 3 '-distal 
conformation (in which the 3' substituent of the outer ring projects down and away 
from the inner ring). The amino-propionic acid and the outer phenolic ring 
assume the transoid conformation, each on opposite sides of the inner ring. The 
torsion angle Xi for the amino- propionic acid is 300*. 

The amino-propionic acid substituent is packed loosely in a polar pocket 
formed by side chains from H2, H4 and S3. The carboxylate group forms direct 
hydrogen bonds with the guanidium group of Arg228 and the amino N of Ser277. 
In addition, Arg262, Arg266 and Asnl79 interact with the carboxylate through 
water-mediated hydrogen bonds. The three arginine residues create a significantly 
positive local electrostatic potential, which may stabilize the negative charge of the 
carboxylate. No hydrogen bond is formed by the amino nitrogen. The 
interactions of the amino-propionic acid substituent are consistent with the fact that 
triac, which lacks the amino nitrogen, has a binding affinity equal to that of T 3 , 
indicating that the amino nitrogen and longer aliphatic chain of T 3 do not 
contribute greatly to binding affinity. 

The diphenyl ether, in contrast, is found buried within the hydrophobic 
core. The inner ring packs in a hydrophobic pocket formed by H3, H5-6, and S3. 
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Pockets for the 3- and 5-methyl' substituents are not completely filled, as expected 
since the van der waals radius of methyl substituent for Dimit is smaller than the 
iodine substituent provided by the thyroid hormone T 3 . Such pockets are typically 
25 to 100 cubic angstroms (although smaller pocket for substitutes are 
5 contemplated in the 40 to 80 cubic angstrom range) and could be filled more 

Lightly with better fitting chemical substitutions, as described herein. 

The outer ring packed tightly in a pocket formed by H3, H5-6, H7, H8, 
Hll and H12, and the loop between H7 and H8. The ether oxygen is found in a 
hydrophobic environment defined by Phe218, Leu287, Leu276, and Leu292. The 

10 absence of a hydrogen bond to the ether oxygen is consistent with its role in 

establishing the correct stereochemistry of the phenyl rings, as suggested by potent 
binding of hormone analogs with structurally similar linkages possessing reduced 
or negligible hydrogen bonding capability. The 3'-isopropyl substituent contacts 
Gly290 and 291. The presence of glycine at this position in the pocket can 

15 explain the observed relationship between activity and the size of S'-substituents. 

Activity is highest for 3'-isopropyl, and decreases with added bulk. The only 
hydrogen bond in the hydrophobic cavity is formed between the phenolic hydroxyl 
and His381 Ne2. The conformation of His381 is stabilized by packing contacts 
provided by Phe405, and Met256. 

20 The presence of a 5' substituent larger than hydrogen affects the binding 

affinity for hormone. The more abundant thyroid hormone, S^^'^'-tetraiodo-L- 
thyronine (T 4 ), contains an iodine at this position, and binds the receptor with 2% 
of the affinity of T 3 . The structure suggests that discrimination against T 4 is 
accomplished through the combination of steric conflict by Met256 and possibly 

25 the constraints imposed by the geometry of the hydrogen bond from His381 to the 

phenolic hydroxyl. The 5' position is a preferred location for introducing a 
chemical modification of C-H at the 5 1 of T3 or and TR agonist, as described 
herein, that produces an extension from the prime ring and results in the creation 
of an antagonist or partial agonist. 

30 Deletion and antibody competition studies suggest the involvement of 

residues Prol62 to Val202 in ligand binding. The region does not directly contact 
hormone in the bound structure, although H2 packs against residues forming the 
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polar pocket that interacts with the amino-propionic acid group. One role for H2, 
then, is to stabilize these residues in the bound state, H2. with 0-strands S3 and 
S4, might also represent a prevalent entry point for ligand, since the amino- 
propionic acid of the ligand is oriented toward this region. Studies of receptor 
binding to T 3 affinity matrices demonstrate that only a linkage to the amino- 
propionic acid is tolerated, suggesting that steric hindrance present in other 
linkages prevent binding. Furthermore, the crystallographic temperature factors 
suggest the coil and 0-strand region is most flexible part of the domain FIG. 7. 
Participation of this region, pan of the hinge domain between the DBD and LBD, 
in binding hormone may provide structural means for ligand binding to influence 
DNA binding, since parts of the Hinge domain contact DNA. 

TR LBD Transcriptional Activation Helix As An Example Of A Nuclear Receptor 
Activation Domain 

In addition to the startling finding that the ligand binding cavity is solvent 
inaccessible when loaded with a ligand, the activation helix of TR LBD presents a 
surface to the ligand cavity for interaction between at least one amino acid and the 
bound ligand. The C-terminal 17 amino acids of the TR, referred to as the 
activation helix or AF-2 (an example of an LBD activation domain), are implicated 
in mediating hormone-dependent transcriptional activation. Although, mutations of 
key residues within the domain decrease ligand-dependent activation it was unclear 
until the present invention whether such mutations directly affected ligand 
coordination. Although some mutations of this domain have been noted to reduce 
or abolish ligand binding, other mutations in more distant sites of the LBD have a 
similar effect. 

Activation domains among nuclear receptors display an analogous three 
dimensional relationship to the binding cavity, which is a region of the LBD that 
binds the molecular recognition domain of a ligand, i.e. the activation domain 
presents a portion of itself to the binding cavity (but necessarily the molecular 
recognition domain of the ligand). Many nuclear receptors are expected to have 
such domains, including the retinoid receptors, RAR and RXR, the glucocorticoid 
receptor GR, and the estrogen receptor ER. Based upon the TR's sequence, the 
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domain is proposed to adopt an amphipathic helical structure. /3-sheet or mixed 
secondary structures, could be present as activation domains in less related nuclear 
receptors. 

Within the activation domain, the highly conserved motif 3><i>XE<M\ where 
5 $ represents a hydrophobic residue, is proposed to mediate interactions between 

the receptors and transcriptional coactivators. Several proteins have been 
identified which bind the TR in a hormone-dependent fashion. One of these, 
Tripl, is related to a putative yeast coactivator Sugl, and also interacts with both 
the C-terminal activation domain and a subset of the basal transcriptional 

10 machinery, suggesting a role in transactivation by the TR. Other proteins, such as 

RIP140, SRC1, (Onate, S.A. et. ah, Science 270:1354-1357 (1995)) and TF-1 (see 
also Ledouarim, B., et. ah. EMBO J. 14:2020-2033 (1995)), also interact with 
other nuclear receptors in a ligand dependent manner through the C-terminal 
domain. Binding of these proteins can be modulated using the TR ligands 

15 described herein especially those TR ligands with extensions that sterically hinder 

the interaction between the highly conserved motif and other proteins. 

The C-terminal activation domain of the TR forms an amphipathic helix, 
H12, which nestles loosely against the receptor to form part of the hormone 
binding cavity. The helix packs with the hydrophobic residues facing inward 

20 towards the hormone binding cavity, and the charged residues, including the 

highly-conserved glutamate, extending into the solvent, as shown in FIG. 8. The 
activation helix of TR LBD presents Phe 401 to the ligand binding cavity and 
permits direct coordination with the hormone i.e. such amino acids interact with 
the ligand forming a van der waals contact with the plane of the outer phenyl ring. 

25 Phe 405 also interacts with His 381, perhaps stabilizing its hydrogen bonding 

conformation, i.e. a favorable hydrogen bond interaction. Participation of Phe 
401 and Phe 405 in binding hormone explains how mutation of these residues 
decreases hormone binding affinity. Furthermore, the impact of these mutations 
on activation likely derives from a role in stabilizing the domain in the bound 

30 structure through increased hydrogen bond interaction of dipole interactions. Glu 

403 extends into the solvent, emphasizing its critical role in transactivation. In its 
observed conformation, presented on the surface as an ordered residue, against a 
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background of predominantly hydrophobic surface, Glu 403 is available to interact 
with activator proteins described herein, as shown in FIG. 9. The other charged 
residues, Glu 405 and Asp 406 are disordered, as the helix frays at Phe 405. 

Two other sequences in the TR, t2 and t3, activate transcription when 
expressed as fusion proteins with a DNA-binding domain. The sequences, 
discovered in the TRB, correspond to TR-or residues Prol58-Ilel68 in HI (r2), 
and Gly290-Leu3 19 in H8 and H9 (t3). Unlike the C-terminal activation domain, 
r2 and t3 do not appear to represent modular structural units in the rat TR-a 
LBD, nor present a surface for protein-protein interactions: the critical 
aspartate/glutamate residues of r3 are located on two separate helices, and do not 
form a single surface; the charged residues of r2 are engaged in ion pair 
interactions with residues of the LBD. Thus, t2 and t3 may not function as 
activation domains in the context of the entire receptor. 

Computational Methods For Designing A Nuclear Receptor LBD Ligand 

The elucidation of the three dimensional structure of a nuclear receptor 
ligand binding domain provides an important and useful approach for designing 
ligands to nuclear receptors using the computational methods described herein. By 
inspecting the FIGURES it can be determined that the nuclear receptor ligand is 
bound in a water inaccessible binding cavity in the LBD and that chemical 
moieties can be added to selected positions on the ligand. Such chemical 
modifications, usually extensions, can fill up the binding cavity represented in the 
FIGURES for a tighter fit (or less water) or can be used to disrupt or make 
contacts with amino acids not in contact with the ligand before the chemical 
modification was introduced or represented in a figure of the three dimensional 
model of the LBD. Ligands that interact with nuclear superfamily members can 
act as agonists, antagonists and partial agonists based on what ligand-induced 
conformational changes take place. 

Agonists induce changes in receptors that place them in an active 
conformation that allows them to influence transcription, either positively or 
negatively. There may be several different ligand-induced changes in the 
receptor's conformation. 
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Antagonists, bind to receptors, but fail to induce conformational changes 
that alter the receptor's transcriptional regulatory properties or physiologically 
relevant conformations. Binding of an antagonist can also block the binding and 
therefore the actions of an agonist. 

Partial agonists bind to receptors and induce only pan of the changes in the 
receptors that are induced by agonists. The differences can be qualitative or 
quantitative. Thus, a partial agonist may induce some of the conformation 
changes induced by agonists, but not others, or it may only induce certain changes 
to a limited extent. 

Ligand-induced Conformational Changes 

As described herein, the unliganded receptor is in a configuration that is 
either inactive, has some activity or has repressor activity. Binding of agonist 
ligands induces conformational changes in the receptor such that the receptor 
becomes more active, either to stimulate or repress the expression of genes. The 
receptors may also have non-genomic actions, some of the known types of 
changes and/or the sequelae of these are listed herein. 
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Heal Shock Protein Binding 

For many of the nuclear receptors ligand binding induces a dissociation of 
heat shock proteins such that the receptors can form dimers in most cases, after 
which the receptors bind to DNA and regulate transcription. 

Nuclear receptors usually have heat shock protein binding domains that 
present a region for binding to the LBD and can be modulated by the b.nding of a 
ligand to the LBD. Consequently, an extended chemical moiety (or more) from 
the ligand that stabilizes the binding or contact of the heat shock protein binding 
domain with the LBD can be designed using the computational methods described 
herein to produce a partial agonist or antagonist. Typically such extended 
chemical moieties will extend past and away from the molecular recognition 
domain on the ligand and usually past the buried binding cavity of the ligand. 

Dimerization and Heterodimerization 

With the receptors that are associated with the hsp in the absence of the 
ligand, dissociation of the hsp results in dimerization of the receptors. 
Dimerization is due to receptor domains in both the DBD and the LBD. Although 
the main stimulus for dimerization is dissociation of the hsp. the ligand-induced 
conformational changes in the receptors may have an additional facilitative 
influence. With the receptors that are not associated with hsp in the absence of 
the ligand, particularly with the TR, ligand binding can affect the pattern of 
dimerization/heterodimerization. The influence depends on the DNA binding site 
context, and may also depend on the promoter context with respect to other 
proteins that may interact with the receptors. A common pattern is to discourage 
monomer formation, with a resulting preference for heterodimer formation over 
dimer formation on DNA. 

Nuclear receptor LBDs usually have dimerization domains that present a 
region for binding to another nuclear receptor and can be modulated by the 
binding of a ligand to the LBD. Consequently, an extended chemical moiety (or 
more) from the ligand that disrupts the binding or contact of the dimerization 
domain can be designed using the computational methods described herein to 
produce a partial agonist or antagonist. Typically such extended chemical moieties 
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will extend past and away from the molecular recognition domain on the ligand 
and usually past the buried binding cavity of the ligand. 



DNA Binding 

5 In nuclear receptors that bind to hsp, the ligand-induced dissociation of hsp 

with consequent dimer formation allows, and therefore, promotes DNA binding. 
With receptors that are not associated (as in the absence of ligand), ligand binding 
tends to stimulate DNA binding of heterodimers and dimers, and to discourage 
monomer binding to DNA. However, with DNA containing only a single half 

10 site, the ligand tends to stimulate the receptor's binding to DNA. The effects are 

modest and depend on the nature of the DNA site and probably on the presence of 
other proteins that may interact with the receptors. Nuclear receptors usually have 
DBD (DNA binding domains) that present a region for binding to DNA and this 
binding can be modulated by the binding of a ligand to the LBD. Consequently, 

15 an extended chemical moiety (or more) from the ligand that disrupts the binding or 

contact of the DBD can be designed using the computational methods described 
herein to produce a partial agonist or antagonist. Typically such extended 
chemical moieties will extend past and away from the molecular recognition 
domain on the ligand and usually past the buried binding cavity of the ligand. 

20 

Repressor Binding 

Receptors that are not associated with hsp in the absence of ligand 
frequently act as transcriptional repressors in the absence of the ligand. This 
appears to be due, in part, to transcriptional repressor proteins that bind to the 
25 LBD of the receptors. Agonist binding induces a dissociation of these proteins 

from the receptors. This relieves the inhibition of transcription and allows the 
transcriptional transactivation functions of the receptors to become manifest. 

Transcriptional Transactivation Functions 
30 Ligand binding induces transcriptional activation functions in two basic 

ways. The first is through dissociation of the hsp from receptors. This 
dissociation, with consequent dimerization of the receptors and their binding to 
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DNA or other proteins in the nuclear chromatin allows transcriptional regulatory 
properties of the receptors to be manifest. This may be especially true of such 
functions on the amino terminus of the receptors. 

The second way is to alter the receptor to interact with other proteins 
involved in transcription. These could be proteins that interact directly or 
indirectly with elements of the proximal promoter or proteins of the proximal 
promoter. Alternatively, the interactions could be through other transcription 
factors that themselves interact directly or indirectly with proteins of the proximal 
promoter. Several different proteins have been described that bind to the receptors 
in a ligand-dependent manner. In addition, it is possible that in some cases, the 
ligand-induced conformational changes do not affect the binding of other proteins 
to the receptor, but do affect their abilities to regulate transcription. 

Nuclear receptors or nuclear receptor LBDs usually have activation 
domains that present a region for binding to DNA and can be modulated by the 
binding of a ligand to the LBD. Consequently, an extended chemical moiety (or 
more) from the ligand that disrupts the binding or contact of the activation domain 
can be designed using the computational methods described herein to produce a 
partial agonist or antagonist. Typically such extended chemical moieties will 
extend past and away from the molecular recognition domain on the ligand and 
usually past the buried binding cavity of the ligand and in the direction of the 
activation domain, which is often a helix as seen in the three dimensional model 
shown in the FIGURES in two dimensions on paper or more conveniently on a 
computer screen. 

Ligand-induced Conformational Change 

Plasma proteins bind hormones without undergoing a conformational 
change through a static binding pocket formed between monomers or domains. 
For example, the tetrameric thyroid-binding plasma protein transthyretin forms a 
solvent-accessible hormone-binding channel at the oligomer interface. The 
structure of the protein is unchanged upon binding hormone with respect to the 
appearance of a buried binding cavity with a ligand bound. 
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However, the structural role for a ligand bound to a nuclear receptor LBD, 
like rat TR-a LBD, predicts that the receptor would differ in the bound and 
unbound states. In the absence of hormone, the receptor would possess a cavity at 
its core, uncharacteristic of a globular protein. A ligand (e.g. hormone) completes 
the hydrophobic core of the active receptor after it binds to the nuclear receptor. 
Ligand binding by the receptor is a dynamic process, which regulates receptor 
function by inducing an altered conformation. 

An exact description of the hormone- induced conformational changes 
requires comparison of the structures of the liganded and the unliganded TR. The 
structure of the unliganded human RXRa may substitute as a model for the 
unliganded TR. The rat TR-a LBD and human RXRa LBDs adopt a similar fold, 
and it is likely that the structural similarity extends to the conformational changes 
after ligand binding. 

There are three major differences between the two structures, which indeed 
appear to be the result of ligand binding. First, the bound rat TR-a LBD structure 
is more compact, with the hormone tightly packed within the hydrophobic core of 
the receptor. By contrast, the unliganded human RXRa LBD contains several 
internal hydrophobic cavities. The presence of such cavities is unusual in folded 
proteins, and is likely a reflection of the unliganded state of the receptor. Two of 
these cavities were proposed as possible binding sites for 9-cis retinoic acid, 
though these multiple sites only partly overlap with the single buried binding 
cavity observed in the liganded rat TR-a LBD. 

The second difference involves Hll in the rat TR-a LBD, which 
contributes part of the hormone binding cavity. Hll, continuous in the rat TR-a 
LBD, is broken at Cys 432 in the RXR, forming a loop between H10 and Hll in 
the hRXRa. This residue corresponds to His381 in the TR, which provides a 
hydrogen bond to the outer ring hydroxy 1 of the ligand. Furthermore, the 
hormone binding cavity occupied by ligand in the rat TR-a LBD is interrupted in 
the hRXRa by the same loop, forming an isolated hydrophobic pocket in the RXR 
with H6 and H7. In the bound rat TR-a LBD, the corresponding helices H7 and 
H8 are contiguous with the binding pocket, and enclose the hormone binding 
cavity from below. 
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The third difference between the two receptors is the position of the C- 
terminal activation domain. While the C-terminal activation domain forms a- 
helices in both receptors, the domain in the rat TR-a LBD follows a proline-rich 
turn, and lies against the receptor to contribute pan of the binding cavity. In 
contrast, the activation domain in the unliganded hRXRa, is pan of a longer helix 
which projects into the solvent. 

These differences lead to a model for an alternate conformation of the TR 
LBD assumed in the absence of ligand. In the unliganded TR, the subdomain of 
the receptor surrounding the hormone binding cavity is loosely packed, with the 
binding cavity occluded by a partly unstructured Hll providing a partial core for 
the receptor. 

Upon binding hormone, residues which form a coil in the unbound receptor 
engage the ligand, and continues Hll. The ordering of Hll could unblock the 
hydrophobic cavity, allowing H7 and H8 to interact with hormone. The extended 
hydrophobic cavity then collapses around the hormone, generating the compact 
bound structure. 

It is possible to predict ligand-induced conformational changes in the C- 
terminal activation domain that rely, in part, on an extended structure in the 
unliganded TR that repacks upon ligand binding. The ligand- induced 
conformation change can be subtle since the amino acid sequence of the rat TR-a 
in the turn (393-PTELFPP-399) significantly reduces the propensity of the peptide 
chain of the rat TR-a to form an a-helix and therefore repacking can be 
accomplished with a minor change in volume. 

After the ligand-induced conformational change occurs, it is likely that the 
conformation of the C-terminal activation domain in the bound structure changes 
packing compared to the unbound form of the receptor. Binding of the ligand 
improves the stability of the activation domain. The activation domain packs 
loosely even in the bound structure, as measured by the distribution of packing 
interactions for the entire LBD. The packing density for the activation domain, 
defined as the number of atoms within 4.5A, is 1.5 standard deviations below the 
mean. For comparison, another surface helix, HI, is 0.5 standard deviations 
below the mean and the most poorly packed pan of the structure, the irregular coil 



WO 97/21993 PCT/US96/20778 

29. 

from residues Ilel96-Asp206, is 2.0 standard deviations below the mean. 
Moreover, the majority of packing contacts for the C-terminal domain in the 
bound receptor are provided either by residues which interact with ligand, such as 
His381, or by the ligand itself. The conformation of these residues can be 
expected to be different in the bound and unbound receptors, and by extension the 
conformation of C-terminal activation domain which relies upon these interactions. 
Without the stabilization provided by a bound ligand, it is likely that the C- 
terminal activation domain is disordered prior to hormone binding. 

The interrelation of ligand-induced conformational changes is evident as 
described herein. For example, His381 from Hll and Phe405 from.H12 interact 
in the bound structure to provide a specific hydrogen bond to the phenolic 
hydroxyl. The ligand-induced changes which affect Hll and H12 are reinforcing, 
and lead to the formation of the compact, bound state. 

Computational Methods Using Three Dimensional Models and 
Extensions Of Ligands 

The three-dimensional structure of the liganded TR receptor is 
unprecedented, and will greatly aid in the development of new nuclear receptor 
synthetic ligands, such as thyroid receptor antagonists. In addition, this receptor 
superfamily is overall well suited to modern methods including three-dimensional 
structure elucidation and combinatorial chemistry such as those disclosed in EP 
335 628, U.S. patent 5,463,564, which are incorporated herein by reference. 
Structure determination using X-ray crystallography is possible because of the 
solubility properties of the receptors. Computer programs that use crystallography 
data when practicing the present invention will enable the rational design of ligand 
to these receptors. Programs such as RASMOL can be used with the atomic 
coordinates from crystals generated by practicing the invention or used to practice 
the invention by generating three dimensional models and/or determining the 
structures involved in ligand binding. Computer programs such as INSIGHT and 
GRASP allow for further manipulation and the ability to introduce new structures. 
In addition, high throughput binding and bioactivity assays can be devised using 
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purified recombinant protein and modern reporter gene transcription assays 
described herein and known in the an in order to refine the activity of a CDL. 

Generally the computational method of designing a nuclear receptor 
synthetic ligand comprises two steps: 

1) determining which amino acid or amino acids of a nuclear receptor LBD 
interacts with a first chemical moiety (at least one) of the ligand using a three 
dimensional model of a crystallized protein comprising a nuclear receptor LBD 
with a bound ligand, and 

2) selecting a chemical modification (at least one) of the first chemical 
moiety to produce a second chemical moiety with a structure to either decrease or 
increase an interaction between the interacting amino acid and the second chemical 
moiety compared to the interaction between the interacting amino acid and the first 
chemical moiety. 

As shown herein, interacting amino acids form contacts with the ligand and 
the center of the atoms of the interacting amino acids are usually 2 to 4 angstroms 
away from the center of the atoms of the ligand. Generally these distances are 
determined by computer as discussed herein and in McRee 1993, however 
distances can be determined manually once the three dimensional model is made. 
Examples of interacting amino acids are described in Appendix 2. See also 
Wagner et al., Nature 378(6558):670-697 (1995) for stereochemical figures of 
three dimensional models. More commonly, the atoms of the ligand and the 
atoms of interacting amino acids are 3 to 4 angstroms apart. The invention can be 
practiced by repeating steps 1 and 2 to refine the fit of the ligand to the LBD and 
to determine a better ligand, such as an agonist. As shown in the FIGURES the 
three dimensional model of TR can be represented in two dimensions to determine 
which amino acids contact the ligand and to select a position on the ligand for 
chemical modification and changing the interaction with a particular amino acid 
compared to that before chemical modification. The chemical modification may 
be made using a computer, manually using a two dimensional representation of the 
three dimensional model or by chemically synthesizing the ligand. The three 
dimensional model may be made using Appendix 2 and the FIGURES. As an 
additional step, the three dimensional model may be made using atomic 
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coordinates of nuclear receptor LBDs from crystallized protein as known in the 
art, see McRee 1993 referenced herein. 

The ligand can also interact with distant amino acids after chemical 
modification of the ligand to create a new ligand. Distant amino acids are 
generally not in contact with the ligand before chemical modification. A chemical 
modification can change the structure of the ligand to make as new ligand that 
interacts with a distant amino acid usually at least 4.5 angstroms away from the 
ligand. Often distant amino acids will not line the surface of the binding cavity 
for the ligand, as they are too far away from the ligand to be part of a pocket or 
surface of the binding cavity. 

The interaction between an atom of a LBD amino acid and an atom of an 
LBD ligand can be made by any force or attraction described in nature. Usually 
the interaction between the atom of the amino acid and the ligand will be the result 
of a hydrogen bonding interaction, charge interaction, hydrophobic interaction, 
van der waals interaction or dipole interaction. In the case of the hydrophobic 
interaction it is recognized that this is not a per se interaction between the amino 
acid and ligand, but rather the usual result, in part, of the repulsion of water or 
other hydrophilic group from a hydrophobic surface. Reducing or enhancing the 
interaction of the LBD and a ligand can be measured by calculating or testing 
binding energies, computationally or using thermodynamic or kinetic methods as 
known in the art. 

Chemical modifications will often enhance or reduce interactions of an 
atom of a LBD amino acid and an atom of an LBD ligand. Steric hinderance will 
be a common means of changing the interaction of the LBD binding cavity with 
the activation domain. Chemical modifications are preferably introduced at C-H, 
C- and C-OH position in ligands, where the carbon is part of the ligand structure 
which remains the same after modification is complete. In the case of C-H, C 
could have 1, 2 or 3 hydrogens, but usually only one hydrogen will be replaced. 
The H or OH are removed after modification is complete and replaced with the 
desired chemical moiety. 

Because the thyroid receptor is a member of the larger superfamily of 
hormone-binding nuclear receptors, the rules for agonist and antagonist 
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development will be recognized by one skilled in the art as useful in designing 
hgands to the entire superfamily. Examining the structures of known agonists and 
antagonists of the estrogen and androgen receptors supports the generality of 
antagonist mechanism of action as shown in FIG. 10. 

The overall folding of the receptor based on a comparison of the reported 
structure of the unliganded RXR and with amino acid sequences of other 
superfamily members reveals that the overall folding of receptors of the 
superfamily is similar. Thus, it is predicted from the structure that there is a 
general pattern of folding of the nuclear receptor around the agonist or antagonist 
-ligand. 

The three dimensional structure of a nuclear receptor with a ligand bound 
leads to the nonobvious observation that a nuclear receptor folds around agonist 
ligands, as the binding cavity fits the agonist, especially the agonist's molecular 
recognition domain, and antagonists commonly have chemical structures that 
extend beyond the ligand, especially the agonist, and would prohibit folding of the 
receptor around the ligand to form a buried binding cavity or other groups that 
have the same effect. The location of the extension could affect the folding in 
various ways as indicated by the structure. Such extensions on antagonists are 
shown in FIG. 10 for various receptors and compared to the corresponding 
agonist. 

For example, an extension towards the carboxy-terminal activation helix 
affects the packing/folding of this helix into the body of the receptor. This in turn 
can affect the ability of this portion of the nuclear receptor to interact with other 
proteins or other portions of the receptor, including transcriptional transactivation 
functions on the opposite end of the linear receptor, or the receptor's amino 
terminus that may interact directly or indirectly with the carboxy-terminal 
transactivation domain (including helix 12). Extensions in this direction can also 
affect the packing of helix 11 of TR (or its analogous helix in nuclear receptors) 
into the body of the receptor and selectively affect dimerization and 
heterodimerization of receptors. An extension pointing towards helix 1 can affect 
the relationship of the DNA binding domain and hinge regions of the receptors 
with the ligand binding domain and selectively or in addition affect the receptors' 
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binding to DNA and/or interactions of receptors with proteins that interact with 
this region of the receptor. Other extensions towards helix 11 can be made to 
affect the packing of this helix and helices 1 and 10 and thereby dimerization. 
Such chemical modifications can be assessed using the computational methods 
described herein. It is also possible that, in some cases, extensions may protrude 
through the receptor that is otherwise completely or incompletely folded around 
the ligand. Such protruding extensions could present a steric blockade to 
interactions with co-activators or other proteins. 

The three dimensional structure with the ligand buried in the binding cavity 
immediately offers a simple description of a nuclear receptor that has a binding 
cavity that contains hinges and a lid, composed of one or more structural 
elements, that move to accommodate and surround the ligand. The ligand to TR 
can be modified on specific sites with specific classes of chemical groups that will 
serve to leave the lid and hinge region in open, partially open or closed states to 
achieve partial agonist or antagonist functions. In these states, the biological 
response of the TR is different and so the structure can be used to design 
particular compounds with desired effects. 

Knowledge of the three-dimensional structure of the TR-T 3 complex leads 
to a general model for agonist and antagonist design. An important novel feature 
of the structural data is the fact that the T 3 ligand is completely buried within the 
central hydrophobic core of the protein. Other ligand-receptor complexes 
belonging to the nuclear receptor superfamily will have a similarly buried ligand 
binding site and therefore this model will be useful for agonist/antagonist design 
for the entire superfamily. 

When design of an antagonist is desired, one needs either to preserve the 
important binding contacts of natural hormone agonist while incorporating an 
"extension group" that interferes with the normal operation of the ligand-receptor 
complex or to generate the requisite binding affinity through the interactions of the 
extensions with receptor domains. 

The model applied to antagonist design and described herein is called the 
"Extension Model." Antagonist compounds for nuclear receptors should contain 
the same or similar groups that facilitate high-affinity binding to the receptor, and 
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in addition, such compounds should contain a side chain which may be large 
and/or polar. This side chain could be an actual extension, giving it bulk, or it 
could be a side group with a charge function that differs from the agonist ligand. 
For example, substitution of a CH 3 for CH : OH at the 21 -position, and alteration at 
the 11 -position from an OH group to a keto group of Cortisol generates 
glucocorticoid antagonist activity (Robsseau, G.G., et. al., J. Mol. Biol. 67:99- 
115 (1972)). However, in most cases effective antagonists have more bulky 
extensions. Thus, the antiglucocorticoid (and antiprogestin) RU486 contains a 
bulky side group at the 11-position (Horwitz, K.B. Endocrine Rev. 13: 146-163 
(1992)). The antagonist compound will then bind within the buried ligand binding 
site of the receptor with reasonably high affinity (100 nM), but the extension 
function will prevent the receptor-ligand complex from adopting the necessary 
conformation needed for transcription factor function. The antagonism (which 
could be in an agonist or antagonist) may manifest itself at the molecular level in a 
number of ways, including by preventing receptor homo/heterodimer formation at 
the HRE, by preventing coactivator binding to receptor monomers, homodimers or 
homo/heterodimers, or by a combination of these effects which otherwise prevent 
transcription of hormone responsive genes mediated by ligand-induced effects on 
the HRE. There are several antagonist compounds for nuclear receptors in the 
prior art (see also Horwitz, K.B., Endocrine Rev. 13:146-163 (1992), Raunnaud 
J.P. et. al., J. Steroid Biochem. 25:811-833 (1986), Keiel S., et. al., Mol. Cell. 
Biol. 14:287-298 (1994) whose antagonist function can be explained by the 
extension hypothesis. These compounds are shown in FIG. 10 along with their 
agonist counterparts. Each of these antagonists contains a large extension group 
attached to an agonist or agonist analogue core structure. Importantly, these 
antagonist compounds were discovered by chance and not designed with a 
structure-function hypothesis such as the extension principle. 

One method of design of a thyroid antagonist using the extension 
hypothesis is provided below as a teaching example. The three-dimensional 
structure of the TR-aDimit complex combined with structure-activity data 
published in the prior art, especially those reference herein, can be used to 
establish the following ligand-receptor interactions which are most critical for 
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high-affinity ligand binding. A physical picture of these interactions is shown in 
FIG. 6. The figure describes the isolated essential contacts for ligand binding. 
Because the ligand is buried in the center of the receptor, the structural spacing 
between these isolated interactions is also important. Thus, our present knowledge 
5 of this system dictates that, for this example, a newly designed ligand for the 

receptor must contain a thyronine structural skeleton, or two substituted aryl 
groups joined by a one-atom spacer. 

The general structure for an antagonist designed by the extension 
hypothesis is exemplified in the following general description of the substituents of 

10 a TR antagonist (referring to Formula 1): RI can have anionic groups such as a 

carboxylate, phosphonate, phosphate, sulfate or sulfite and is connected to the ring 
with a 0 to 3 atom linker, comprising one or more C, O, N, S atoms, and 
preferably a 2 carbon linker. Such Rl can be optionally substituted with an amine 
(e.g. -NH2). R3 and R5 are small hydrophobic groups such as -Br, -I, or -CH3. 

15 R3 and R5 can be the same substituents or different. R 3 ' can be a hydrophobic 

group that may be larger than those of R3 and R5, such as -I, -CH3, -isopropyl, - 
phenyl, -benzyl, 5 and 6 ring heterocycles. R/ is a group that can participate in a 
hydrogen bond as either a donor or acceptor. Such groups include -OH, -NH 2 , 
and -SH. R 5 ' is an important extension group that makes this compound an 

20 antagonist. R 5 ' can be a long chain alkyl (e.g. 1 to 9 carbons, straight chain or 

branched), aryl (benzyl, phenyl and substituted benzyl and phenyl rings (e.g with 
halogen, alkyl (1 and 5 carbons) and optionally connected to the ring by a -CH2-), 
heterocycle (e.g. 5 or 6 atoms, preferably 5 carbons and 1 nitrogen, or five 
carbons), which can optionally include polar (e.g. -OH, -NH 2 , and -SH), cationic 

25 (e.g. -NH3, N(CH)3), or anionic (carboxylate, phosphonate, phosphate or sulfate) 

groups. R 5 ' can also be a polar (e.g. -OH, -NH 2 , and -SH), cationic (e.g. -NH3, - 
N(CH3)3), and anionic (carboxylate, phosphonate, phosphate or sulfate) groups. 
X is the spacer group that appropriately positions the two aromatic rings. This 
group is usually a one-atom spacer, such as O, S, SO, S02, NH, NZ where Z is 

30 an alkyl, CH2, CHOH, CO, C(CH3)OH, and C(CH3)(CH3). R2, R6, R2' and 

R6' can be -F, and -CI and are preferably H. 
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A TR ligand can also be described as a substituted phenylated 3,5 diiodo 
tyrosine with substituted R5' and R3' groups. R5' can be a long chain alkyl (e.g. 
4 to 9 carbons, straight chain or branched), aryl (benzyl, phenyl and substituted 
benzyl and phenyl rings (e.g with halogen, alkyl (1 and 5 carbons) and optionally 
connected to the ring by a -CH2-), heterocycle (e.g. 5 or 6 atoms, preferably 5 
carbons and 1 nitrogen, or five carbons), which can optionally include polar (e.g. 
-OH, -NH 2 , and -SH), cationic (e.g. -NH3, N(CH)3), or anionic (carboxylate, 
phosphonate, phosphate or sulfate) groups. R5' can also be a polar (e.g. -OH, 
-NH 2 , and -SH), cationic (e.g. -NH3, N(CH)3), and anionic (carboxylate, 
phosphonate, phosphate or sulfate) groups. R3' can be -IsoPr, halogen, -CH3, 
alkyl (1 to 6 carbons) or aryl (benzyl, phenyl and substituted benzyl and phenyl 
rings (e.g with halogen, alkyl (1 and 5 carbons) and optionally connected to the 
ring by a -CH2-), heterocycle (e.g. 5 or 6 atoms, preferably 5 carbons and 1 
nitrogen, or five carbons), which can optionally include polar (e.g. -OH, -NH>, 
and -SH), cationic (e.g. -NH3, N(CH)3), or anionic (carboxylate, phosphonate, 
phosphate or sulfate) groups. 

A TR antagonist can also be a modified T 3 agonist (having a diphenyl 
structure) wherein R 5 ' is alkyl, aryl, 5- or 6-membered heterocyclic aromatic, 
heteroalkyl, heteroaryl, arylalkyl, heteroaryl alkyl, polyaromatic, 
polyheteroaromatic, polar or charged groups, wherein said R s ' may be substituted 
with polar or charged groups. The R5' groups are defined, as described herein. 

Using these methods the ligands of this example preferably have the 
following properties: 

1 . The compounds should bind to the TR with high affinity (for example 
100 nM). 

2. The compounds should bind the receptor in the same basic orientation 
as the natural hormone. 

3. The extension group R5' should project toward the activation helix 
(C-terminal helix) of the receptor. 

4. The appropriate substituent at R5' should perturb the activation helix 
from its optimal local structure needed for mediating transcription. 
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Antagonists may also be designed with multiple extensions in order to 
block more than one aspect of the folding at any time. 

TR ligands (e.g. super agonists) can be designed (and synthesized) to 
enhance the interaction of at least one amino acid with at least one chemical 

5 moiety on the ligand ? s molecular recognition domain. One method is to enhance 

the charge and polar interactions by replacing the carboxylate of T 3 (Rl position) 
with phosphonate, phosphate, sulfate or sulfite. This enhances the interaction with 
Arg 262, Arg 266 and Arg 228. The interaction of at least one amino acid with at 
least one chemical moiety on the ligand's molecular recognition domain can also 

10 be enhanced by increasing the size of Ri group to fill the space occupied by water 

when Dimit is bound (referring to Rl). Preferably the group has a 
complementary charge and hydrophobicity to the binding cavity. 

Another way of improving the interaction of at least one amino acid with at 
least one chemical moiety on the ligand's molecular recognition domain is to 

15 restrict the conformation of the dihedral angle between the two phenyl rings of the 

thyronine ligand in solution. In solution the planes of two phenyl rings are 
orthogonal where the dihedral angle is 90°. In the TR Dimit structure, the 
dihedral angle is close to 60°. A TR ligand design that fixes the angle between 
the two phenyl rings will lead to tighter binding. Such a ligand may be made by 

20 connecting the R6' and the R5 positions of a thyronine or a substituted thyronine- 

like diphenyl. The size of the cyclic connection can fix the angle between the two 
phenyl rings. Referring specifically to Formula 1, the following cyclic 
modifications are preferred: 1) 1*5 is connected to R^\ 2) R 3 is connected to R 2 ' or 
3) R 5 is connected to R^ and R 3 is connected to R2\ The connections can be 

25 made by an alkyl or heteroalkyl chain having between 1 to 6 atoms and preferably 

from 2 to 4 carbon atoms or other atoms. Any position of the heteroalkyl chain 
can be N, O, P or S. The S and P heteroatoms along said heteroalkyl chain are in 
any of their possible oxidative states. The N heteroatom or any carbon along the 
alkyl or heteroalkyl chain may have one or more Z substituents, wherein Z is 

30 alkyl, heteroalkyl, aryl, heteroaryl, 5- or 6-membered heterocyclic aromatic. 

These compounds can be claimed with the proviso that Formula 1 does not include 
any prior art compound as of the filing date of this application. 
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The interaction of at least one amino acid with at least one chemical moiety 
on the ligand's molecular recognition domain can also be enhanced by selecting a 
chemical modification that fills the unfilled space between a TR ligand and the 
LBD in the area of the bridging oxygen (such as in T3, triac or Dimit). Thus, a 
slighter larger moiety that replaces the ether oxygen can enhance binding. Such a 
linker may be a mono- or geminal- disubstituted carbon group. A group 
approximately the same size as oxygen but with greater hydrophobicity is 
preferred as well as small, hydrophobic groups for the disubstituted carbon. 

TR-ar AND TR-0 SELECTIVITY FOR THE THYROID HORMONE RECEPTOR 

Using the method described herein ligands can be designed that selectively 
bind to the alpha more than the beta TR. The X-ray crystallographic structure of 
the rat TR-a LBD provides insight into design of such ligands. 

The three dimensional structure reveals that the major difference between 
the TR-a and TR-0 in the ligand binding cavity resides in amino acid Ser 277 
(with the side group -CH20H) in the rat TR-a and whose corresponding residue is 
331, asparagine (with the side group -CH2CONH2), in the human TR-jS. The side 
chain in human TR-0 is larger, charged and has a different hydrogen bonding 
potential, which would allow the synthesis of compounds that discriminate between 
this difference. 

For example, in the complex of TRa with triac, Ser277 does not participate 
in ligand binding. The absence of a role for Ser277 (Asn331 in beta) is consistent 
with the equal affinity of triac for the alpha and beta isoforms, and indirectly 
supports the contention that alpha/beta selectivity resides in the amino acid 
substitution Ser277 to Asn331 and its interaction with Arg228. 

In terms of ligand design, these differences mean that for 0-selective 
ligands, some or all of the following differences should be exploited: 

1. The presence of a larger side chain asparagine. 

2. The ability of the carbonyl group on the side chain to provide a 
strong hydrogen bond acceptor. 

3. The ability of the amido group on the side chain to provide a two 
hydrogen bond donors. 
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4. Adjustment of polarity to reorganize the trapped water in the T3 
pocket. 

In terms of pharmaceutical design, these differences mean that for 
a-selective ligands, some or all of the following differences should be exploited: 
5 1. The presence of a smaller side group. 

2. The ability of the hydroxy 1 on the -CH20H side group carbonyl 
group on the side chain to provide a weak hydrogen donor. 

3. Adjustment of polarity to reorganize the trapped water in the T3 
pocket. 

10 In both cases these differences can be exploited in a number of ways. For 

example, they can also be used with a software set for construction of novel 
organic molecules such as LUDI from Biosym-MSI. 

Methods of Treatment 
15 The compounds of Formula 1 can be useful in medical treatments and 

exhibit biological activity which can be demonstrated in the following tests: 

(i) the induction of mitochondrial a-glycerophosphate dehydrogenase 

(GPDH:EC 1.1.99.5). This assay is particularly useful since in certain species 

e.g. rats it is induced specifically by thyroid hormones and thyromimetics in a 
20 close-related manner in responsive tissues e.g. liver, kidney and the heart 

(Westerfield, W.W., Richert, D.A. and Ruegamer, W.R., Endocrinology, 1965, 

77, 802). The assay allows direct measurement in rates of a thyroid hormone-like 

effect of compounds and in particular allows measurement of the direct thyroid 

hormone-iike effect on the heart; 
25 (ii) the elevation of basal metabolic rate as measured by the increase in 

whole body oxygen consumption (see e.g., Barker et al., Ann. N. Y. Acad. Sci.. 

86:545-562 (I960)); 

(iii) the stimulation of the rate of beating of atria isolated from animals 

previously dosed with thyromimetrics (see e.g., Stephan et al., Biochem. 
30 Pharmacol. (1992) 13:1969-1974; Yokoyama et al., J. Med. Chem. 38:695-707 

(1995)); 
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(iv) the change in total plasma cholesterol levels as determined using a 
cholesterol oxidase kit (for example, the Merck CHOD iodine colorimetric kit. 
see also, Stephan et al. (1992)); 

(v) the measurement of LDL (low density lipoprotein) and HDL (high 
density lipoprotein) cholesterol in lipoprotein fractions separated by 
ultracentrirugation; and p (vi) the change in total plasma triglyceride levels as 
determined using enzymatic color tests, for example the Merck System GPO-PAP 
method. 

The compounds of Formula 1 can be found to exhibit selective 
thyromimetic activity in these tests, 

(a) by increasing the metabolic rate of test animals, and raising hepatic 
GPDH levels at doses which do not significantly modify cardiac GPDH levels. 

(b) by lowering plasma cholesterol and triglyceride levels, and the ratio 
of LDL to HDL cholesterol at doses which do not significantly modify cardiac 
GPDH levels. 

The compounds of Formula 1 may therefore be used in therapy, in the 
treatment of conditions which can be alleviated by compounds which selectively 
mimic the effects of thyroid hormones in certain tissues whilst having little or no 
direct thyromimetic effect on the heart. For example, compounds of Formula 1 
which raise hepatic GPDH levels and metabolic rate at doses which do not 
significantly modify cardiac GPDH levels are indicated in the treatment of obesity. 

Agonists of Formula 1 will lower total plasma cholesterol, the ratio of 
LDL-cholesterol to HDL-cholesterol and triglyceride levels at doses which do not 
significantly modify cardiac GPDH levels are indicated for use as general 
antihyperlipidaemic (antihyperlipoproteinaemic) agents i.e. in the treatment of 
patients having elevated plasma lipid (cholesterol and triglyceride) levels. In 
addition, in view of this effect on plasma cholesterol and triglyceride, they are 
also indicated for use as specific anti-hypercholesterolemic and anti- 
hypertriglyceridaemic agents. 

Patients having elevated plasma lipid levels are considered at risk of 
developing coronary heart disease or other manifestations of atherosclerosis as a 
result of their high plasma cholesterol and/or triglyceride concentrations. Further, 
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since LDL-cholesterol is believed to be the lipoprotein which induces 
atherosclerosis, and HDL-cholesterol believed to transport cholesterol from blood 
vessel walls to the liver and to prevent the build up of atherosclerotic plaque, anti- 
hyperlipidemic agents which lower the ratio of LDL-cholesterol to HDL 
5 cholesterol are indicated as anti-atherosclerotic agents, herein incorporated by 

reference U.S. patents 4,826,876 and 5,466,861. 

The present invention also provides a method of producing selective 
thyromimetic activity in certain tissues except the heart which comprises 
administering to an animal in need thereof an effective amount to produce said 
10 activity of a compound of Formula 1 or a pharmaceutical^ acceptable salt thereof. 

The present invention also relates to a method of lowering plasma lipid 
levels and a method of lowering the ratio of LDL-cholesterol to HDL-cholesterol 
levels by suitably administering a compound of this invention or a 
pharmaceutically acceptable sale thereof. 
15 In addition, compounds of Formula 1 may be indicated in thyroid hormone 

replacement therapy in patients with compromised cardiac function. 

In therapeutic use the compounds of the present invention are usually 
administered in a standard pharmaceutical composition. 

The present invention therefore provides in a further aspect pharmaceutical 
20 compositions comprising a compound of Formula 1 or a pharmaceutically 

acceptable salt thereof and a pharmaceutically acceptable carrier. Such 
compositions include those suitable for oral, parenteral or rectal administration. 



Pharmaceutical Compositions 
25 Compounds of Formula 1 and their pharmaceutically acceptable salts which 

are active when given orally can be formulated as liquids for example syrups, 

suspensions or emulsions, tablets, capsules and lozenges. 

A liquid composition will generally consist of a suspension or solution of 

the compound or pharmaceutically acceptable salt in a suitable liquid carrier(s), for 
30 example ethanol, glycerine, sorbitol, non-aqueous solvent such as polyethylene 

glycol, oils or water, with a suspending agent, preservative, surfactant, wetting 
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agent, flavoring or coloring agent. Alternatively, a liquid formulation can be 
prepared from a reconstitutable powder. 

For example a powder containing active compound, suspending agent, 
sucrose and a sweetener can be reconstituted with water to form a suspension; and 
a syrup can be prepared from a powder containing active ingredient, sucrose and a 
sweetener. 

A composition in the form of a tablet can be prepared using any suitable 
pharmaceutical carrier(s) routinely used for preparing solid compositions. 
Examples of such carriers include magnesium stearate, starch, lactose, sucrose, 
microcrystalline cellulose and binders, for example polyvinylpyrrolidone. The 
tablet can also be provided with a color film coating, or color included as part of 
the carrier(s). In addition, active compound can be formulated in a controlled 
release dosage form as a tablet comprising a hydrophilic or hydrophobic matrix. 

A composition in the form of a capsule can be prepared using routine 
encapsulation procedures, for example by incorporation of active compound and 
excipients into a hard gelatin capsule. Alternatively, a semi-solid matrix of active 
compound and high molecular weight polyethylene glycol can be prepared and 
filled into a hard gelatin capsule; or a solution of active compound in polyethylene 
glycol or a suspension in edible oil, for example liquid paraffin or fractionated 
coconut oil can be prepared and filled into a soft gelatin capsule. Compound of 
Formula 1 and their pharmaceutical ly acceptable salts which are active when given 
parenterally can be formulated for intramuscular or intravenous administration. 

A typical composition for intra-muscular administration will consist of a 
suspension or solution of active ingredient in an oil, for example arachis oil or 
sesame oil. A typical composition for intravenous administration will consist of a 
sterile isotonic aqueous solution containing, for example active ingredient, 
dextrose, sodium chloride, a co-solvent, for example polyethylene glycol and, 
optionally, a chelating agent, for example ethylenediamine tetracetic acid and an 
anti-oxidant, for example, sodium metabisulphite. Alternatively, the solution can 
be freeze dried and then reconstituted with a suitable solvent just prior to 
administration. 
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Compounds of structure'(l) and their pharmaceutical ly acceptable salts 
which are active on rectal administration can be formulated as suppositories. A 
typical suppository formulation will generally consist of active ingredient with a 
binding and/or lubricating agent such as a gelatin or cocoa butter or other low 
5 melting vegetable or synthetic wax or fat. 

Compounds of Formula 1 and their pharmaceutical ly acceptable salts which 
are active on topical administration can be formulated as transdermal 
compositions. Such compositions include, for example, a backing, active 
compound reservoir, a control membrane, liner and contact adhesive. 
10 The typical daily dose of a compound of Formula 1 varies according to 

individual needs, the condition to be treated and with the route of administration. 
Suitable doses are in the general range of from 0.001 to 10 mg/kg bodyweight of 
the recipient per day. 

Within this general dosage range, doses can be chosen at which the 
15 compounds of Formula 1 lower plasma cholesterol levels and raise metabolic rate 

with little or no direct effect on the heart. In general, but not exclusively, such 
doses will be in the range of from 0.5 to 10 mg/kg. 

In addition, within the general dose range, doses can be chosen at which 
the compounds of Formula 1 lower plasma cholesterol levels and have little or no 
20 effect on the heart without raising metabolic rate. In general, but not exclusively, 

such doses will be in the range of from 0.001 to 0.5 mg/kg. 

It is to be understood that the 2 sub ranges noted above are not mutually 
exclusive and that the particular activity encountered at a particular dose will 
depend on the nature of the compound of Formula 1 used. 
25 Preferably, the compound of Formula 1 is in unit dosage form, for 

example, a tablet or a capsule so that the patient may self-administer a single dose. 
In general, unit doses contain in the range of from 0.05-100 mg of a compound of 
Formula 1. Preferred unit doses contain from 0.05 to 10 mg of a compound of 
Formula 1. 

30 The active ingredient may be administered from 1 to 6 times a day. Thus 

daily doses are in general in the range of from 0.05 to 600 mg per day. 
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Preferably, daily doses are in the range of from 0.05 to 100 mg per day. Most 
preferably from 0.05 to 5 mg per day. 

EXAMPLES 
Example 1 - Synthesis of TR Ligands 

Many TR ligands are known in the art, including T4 (thyroxine), T3, T2 
and TS-9. See Jorgensen, Thyroid Hormones and Analogs, in 6 Hormonal 
Proteins and Peptides, Thyroid Hormones 107-204 (Choh Hao Li ed., 1978), 
incorporated by reference herein.. 

The syntheses of several TR ligands are described below. 

Synthesis of TS1, TS2, TS3, TS4, TS5 

TS1, TS2, TS3, TS4 and TS5 and analogs thereof can all be prepared by 
simple acylation of the nitrogen atom of any thyronine analog, including T3 
(3,5,3'-triiodo-L-thyronine), T4 (thyroxine) and 3,5-diiodothyronine. TS1 and 
TS2 are synthesized by reacting T3 with Ph 2 CHC0 2 NHS (N-hydroxy succinimide- 
2,2-diphenylacetate) and C 16 H 33 C0 2 NHS, respectively. TS3 is synthesized by 
reacting T3 with FMOC-C1 (fluorenylmethyloxycarbonylchloride). TS4 is 
synthesized by reacting T3 with tBGC 2 0 (tBOC anhydride or di-t- 
butyldicarbonate). TS5, which differs from TS1-4 by having a -H instead of an -I 
at the R* 3 position, is synthesized by reacting 3,5-diiodothyronine with tBOC 2 0. 
The general reaction scheme for TS1, TS2, TS3, TS4 and TS5 is depicted in FIG, 
11. It should be noted that in the reaction scheme, both TS5 and its precursor 
both have a hydrogen rather than an iodine at the R* 3 position. 

Synthesis of TS6 and TS7 

TS6 is synthesized by reacting TS5 with paranitrophenylisocyanate. TS7 is 
synthesized by reacting TS6 with TFA (trifluoroacetic acid), which cleaves the 
tBOC group. These reactions are simple organic synthesis reactions that can be 
performed by anyone of ordinary skill in the art. The synthetic scheme for TS6 
and TS7 is diagrammed in FIG. 12. 
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Synthesis of TS8 

TS8 is synthesized by reacting TS5 with Ph,CHNH 2 (diphenylmethylamine) 
in the presence of triethylamine and any amide forming condensing reagent, such 
as TBTU (hydroxybenztriazoleuronium tetrafluoroborate) or HBTU 
(hydroxybenztriazoleuronium hexafluorophosphate). The synthesis scheme for 
TS8 is depicted in FIG. 13. 



Synthesis of 3,5-dhodo-3'isopropylthyronine derivatives 

For designing a class of antagonists, it is important to have a hydrophobic 

10 group at the 3' position as well as an extension at the 5' position. Preferred 

hydrophobic groups at the 3' position include: methyl, benzyl, phenyl, iodo, and 
heterocyclic structures. The synthesis of a S.S-diiodo-S'-isopropyl-S'-substituted 
thyronine is described below. The example provided describes the specific steps 
for synthesizing the TS10 compound, but this general reaction scheme can be used 

15 by one of ordinary skill in the art to synthesize any number of 3,5,-diiodo-3'- 

isopropy 1-5 '-substituted thyronine derivatives, which are characterized by having 
an extension at the 5' position. Additional compounds of this class can be 
synthesized using known organic synthesis techniques. 

The synthesis of TS10 is described below and is depicted in FIG. 14. 

20 Numbers used in the reaction scheme for TS10 indicating the reaction product for 

each step are in parentheses. 

2-Formvl-6- i<sn prnpvlanisole (1) : 2-formyl-6-isopropylanisole (10.0 g, 61 
mmol), as made by Casiraghi, et al. JCS Perkin I, 1862 (1980) (incorporated by 
reference), is added dropwise to a suspension of sodium hydride (3.7 g, 153 

25 mmol) in 50 mL THF and 50 mL of DMF in a round bottom flask. The addition 

generates an exothermic reaction and formation of a gray solid. Methyl iodide 
(26.0 g, 183 mmol) is then added dropwise and the reaction mixture is stirred at 
room temperature for 5 hours. The reaction mixture is quenched with 20 mL of 
water, then poured into 500 mL of water, and is extracted with ether (2 x 300 

30 mL). The ether layers are combined, washed with water (5 x 1000 mL), dried 

over magnesium sulfate and concentrated in vacuo to provide 10.2 g (94%) of the 
title compound, with the following 'H NMR (CDC1,) properties: d 10.30 (s, 1H), 
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7.63 (d, 1H, J=3 Hz). 7.50 (d, IH, J=3 Hz), 7.13 (t, 1H. 1=3 Hz), 3.81 (s, 
3H), 3.31 (heptet. IH. J = 7.5 Hz), 1.19 (d, 6H, J=7.5 Hz). 

2-(2-Hvdroxvnonvl)-6- isoproDvlani.sole (not shown in scheme)- 
Octylmagnesium chloride (8.4 mL, 16.9 mmol, 2.0 M) is added dropwise to a 
solution of 1 (1.5 g, 8.4 mmol) in 10 mL THF at -78° C. The reaction mixture is 
stirred for 2 hours with warming to room temperature. The reaction mixture is 
diluted with 50 mL ether and poured into 50 mL water. The ether layer is washed 
with brine (1 x 50 mL), dried over sodium sulfate, and concentrated in vacuo. 
Flash chromatography (silica gel, 10% ether/hexane - 15% ether/hexane) 
provides 734 mg (30%) of the title compound with the following 'H NMR (CDC1,) 
properties: d 7.33-7.10 (m, 3H), 5.00 (br. s, 1H), 3.81 (s, 3H), 3.33 (heptet, 
IH, J = 7 Hz) 1.90-1.19 (m, 14H), 0.86 (t, 3H, J=6.5 Hz); HRMS (EI), found: 
292.2404; calc'd: 292.2402. 

2-nonyl-6-isopropvlanisole (2)- Compound 2 (663 mg, 2.3 mmol) is 
dissolved in solution of 5 mL ethanol and 5 mL acetic acid, and a spatula tip of 
palladium on carbon catalyst is added. The reaction mixture is then charged with 
hydrogen gas (using a simple balloon and needle) and the mixture is stirred at 
room temperature overnight. The next day, the reaction mixture is poured into 
ether (100 mL) and the ether layer is extracted with saturated sodium bicarbonate 
(3 x 100 mL). The ether layer is dried over sodium sulfate and concentrated in 
vacuo to provide 581 mg (91 %) of (2) with the following 'H NMR (CDC1 3 ) 
properties: d 7.14-7.00 (m, 3H), 3.75 (s. 3H), 3.36 (heptet, IH, J=6.8 Hz), 2.63 
(t, 2H, J=7.5 Hz), 1.68-1.15 (m, 14H), 0.86 (t, 3H, J=5.5 Hz); HRMS (EI), 
mass found: 276.2459; calculated: 276.2453. 

Thyronine adduct (4): Fuming nitric acid (0.071 mL) is added to 0.184 
mL acetic anhydride chilled to -5° C. Iodine (66 mg) is added to this mixture 
followed by trifluoroacetic acid (0.124 mL). This mixture is stirred for 1 hour 
with warming to room temperature, at which point all of the iodine is dissolved. 
The reaction mixture was then concentrated in vacuo to provide an oily semi-solid 
material. The residue was dissolved in 0.7 mL of acetic anhydride and cooled to - 
20° C. A solution of anisole (2) (581 mg. 2.1 mmol) in 1.2 mL acetic anhydride 



WO 97/2,993 PCT/US96/20778 

47. 

and 0.58 mL TFA is added dropwise. The reaction mixture is stirred at -20° for 
1 hour, then stirred overnight with warming to room temperature. The reaction 
mixture is partitioned between water and methylene chloride. The methylene 
chloride layer is dried over sodium sulfate and concentrated in vacuo to provide 

5 the iodonium salt (3) as an oil. This material is not purified or characterized, and 

is directly introduced into the coupling reaction. 

N-Trifluoroacetyl-3,5-diiodotyrosine methyl ester (552 mg, 1.0 mmol) 
prepared according to the procedure of N. Lewis and P. Wallbank, Synthesis 1103 
(1987) (incorporated by reference) and all of the crude iodonium salt (3) from 

10 above is dissolved in 5 mL of anhydrous methanol. Diazabicyclo[5.4.0]undecane 

(DBU) (183 mg, 1.2 mmol) and a spatula tip of copper-bronze are added and the 
resulting mixture is stirred at room temperature overnight. The next day, the 
reaction mixture is filtered, and the filtrate is concentrated in vacuo. The crude 
residue is purified by flash chromatography (silica gel, 10% ethyl acetate/hexane) 

15 to provide 30 mg (4%) of the protected thyronine adduct (4). 

Deprotected thyronine (TS10) : The protected thyronine 4 (30 mg, 0.04 
mmol) is dissolved in a mixture of 2.25 mL acetic acid and 2.25 mL 49% 
hydrobromic acid. The reaction mixture is heated to reflux for 5 hours. The 
reaction mixture is cooled to room temperature, and the solvents are removed in 

20 vacuo. Water is added to triturate the oily residue into a gray solid. This solid 

material is filtered, washed with water, and dried over P,O s in vacuo to provide 24 
mg (81%) of the title compound, TS10, with the following 'H NMR (CDCl 3 ) 
properties: d 7.57 (s, 1H), 6.86 (s, 1H), 6.45 (s, 1H), 6.34 (s, 1H), 4.81 (m, 
1H), 3.86 (s, 3H), 3.71 (s, 3H), 3.33-3.05 (m, 3H), 2.58-2.47 (m, 2H), 1.62-0.76 

25 (m, 23H); MS (LSIMS): M + =.817.0. 

As mentioned above, this reaction scheme can be modified by one of 
ordinary skill in the art to synthesize a class of compounds characterized by 3,5- 
diiodo-3'isopropylthyronine derivatives, wherein (1) the 3* isopropyl group can be 
replaced with a hydrophobic group, including methyl, benzyl, phenyl, iodo, and 

30 heterocyclic structures, and (2) a wide variety of chemical structures can be 

incorporated at the 5' position, including alkyl groups, planar aryl, heterocyclic 
groups, or polar and/or charged groups. 
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The aldehyde (1) in the above reaction scheme is a versatile synthetic 



intermediate which allows for the attachment of a variety of chemical moieties to 
the 5' position of the final thyronine derivative. In addition, a variety of chemical 
reactions can be used to attach the chemical moieties. These reactions are well 
known in the an and include organometallic additions to the aldehyde (including 
Grignard reagents, organolithiums, etc.), reductive amination reactions of the 
aldehyde with a primary or secondary amine, and Wittig olefination reactions with 
a phosphorous ylid or stabilized phosphonate anion. Other possibilities include 
reduction of the aldehyde to a benzyl alcohol allowing for etherification reactions 
at the 5" position. As mentioned above, these methods allow for a wide variety of 
chemical structures to be incorporated at the 5' position of the final thyronine 
derivative, including alky! groups, planar aryl, heterocyclic groups or polar and/or 
charged groups. 

Synthesis of 3, 5-dibromo-4-(3\5'-diisopropyl-4'-hydroxyphenoxy) benzoic acid 
(Compound 11). 



(a) A mixture of 2,6-diisopropyl phenol (20 g, 0.11 mol), potassium 
carbonate (62 g, 0.45 mol), acetone (160 ml) and methyl iodide (28 ml, 0.45 
mole) is refluxed for three days. The reaction mixture is filtered through celite, 
evaporated, dissolved in ether, washed twice with 1M sodium hydroxide, dried 
over magnesium sulphate and concentrated to afford 15.1 g (0.08 mol, 70%) of 
2,6-diisopropyl anisole as a slightly yellow oil. 

(b) Fuming nitric acid (12.4 ml, 265 mmol) is added dropwise to 31.4 ml 
of acetic anhydride which is cooled in a dry ice/carbon tetrachloride bath. Iodine 
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11.3 g, 44.4 mmol) is added in one portion followed by dropwise addition of 
trifluoroacetic acid (20.5 ml, 266 mmole). The reaction mixture is stirred at room 
temperature until all the iodine is dissolved. Nitrogen oxides are removed by 
flushing nitrogen into the vessel. The reaction mixture is concentrated, the residue 
is dissolved in 126 ml of acetic anhydride and is cooled in a dry ice/carbon 
tetrachloride bath. To the stirred solution 2,6-diisopropylanisole (51 g, 266 mmol) 
in 150 ml of acetic anhydride and 22.6 ml of trifluoroacetic acid is added 
dropwise. The reaction mixture is left to stand at room temperature over night 
and then is concentrated. The residue is taken up in 150 ml of methanol and 
treated with 150 ml of 10% aqueous sodium bisulfite solution and 1 liter of 2M 
sodium borotetrafluoride solution. After the precipitate aggregates, petroleum 
ether is added and the supernatant is decanted. The precipitate is triturated with 
petroleum ether, filtered, washed with petroleum ether and dried at room 
temperature in vacuo. This affords 34 g (57 mmol, 65%) of bis(3,5-diisopropyl-4- 
methoxyphenyl)iodonium tetrafluoroborate as a white solid. 

(c) To a stirred solution of 3,5-dibromo-4-hydroxybenzoic acid (12 g, 40.5 
mmol) in 250 ml of methanol, thionyl chloride (3 ml) is added dropwise. The 
reaction mixture is refluxed for five days, water is added and the precipitated 
product is filtered off. The residue is dissolved in ethyl acetate. From the 
aqueous phase, methanol is removed by concentration. The aqueous phase is then 
saturated with sodium chloride, and extracted with ethyl acetate. The combined 
organic phases are dried over magnesium sulphate, filtered and concentrated. This 
gives 12.5 g (40.5 mmol, 100%) of 3,5-dibromo-4-hydroxymethyl benzoate as a 
white crystalline solid. 

(d) The products obtained in steps b and c are reacted with each other 
according to the following protocol. To bis(3,5-diisopropyl-4- 
methoxyphenyl)iodonium tetrafluoroborate (2.86 g, 4.8 mmole) and copper bronze 
(0.42 g, 6.4 mmole) in 7 ml. of dichloromethane at 0° C is added dropwise a 
solution of 3,5-dibromo-4-hydroxymethyl benzoate (1.0 g, 3.2 mmole) and 
triethylamine (0.36 g, 3.5 mmole) in 5 ml of dichloromethane. The reaction 
mixture is stirred in the dark for eight days and then is filtered through celite. 
The filtrate is concentrated and the residue is purified by column chromatography 
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(silica gel, 97:3 petroleum ether/ethyl acetate) to give 0.62 g (1.2 mmole, 39%) of 
3,5-dibromo-4-(3' ( 5'-diisopropyl-4'-methoxyphenoxy)methylbenzoate as a solid. 

(e) The product from step d (0.2 g, 0.4 mmole) is dissolved in 2 ml. 
dichloromethane, is put under nitrogen and is cooled at -40° C. To the stirred 
solution is added 1M BBr 3 (1.2 ml, 1.2 mmole) dropwise. The reaction mixture is 
allowed to reach room temperature and then is left over night. It is cooled to 0° 
C and then hydrolyzed with water. Dichloromethane is removed by concentration 
and the aqueous phase is extracted with ethyl acetate. The organic phase is 
washed with 1M hydrochloric acid and brine. Then it is dried over magnesium 
sulphate, filtered and concentrated. The residue is chromatographed (silica, 
96:3.6:0.4 dichloromethane/methanol/acetic acid) producing 93 mg (0.2 mmole, 
51%) of 3,5-dibromo-4-(3',5'-diisopropyl-4'-hydroxyphenoxy)benzoicacid as a 
white solid. 'H nmr (CDC1 3 ) 8 1.23 (d, 12H, methyl), 3.11 (m, 2H, CH), 6.50 
(s, 2H, 2,6-H) 8.33 (s, 2H, 2',6 , -H). 

TABLE 1 and FIG. 15 depict the structures of several TR ligands. 
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Example 2 - Receptor Binding Assays of TR Ligands 

To test the ability of synthesized TR ligands to bind to a thyroid receptor (TR), 
the binding affinity of a TR ligand for TR is assayed using TR's prepared from rat liver 
nuclei and 125, T, as described in J.D. Apriletti, J.B. Baxter, and T.N. Lavin, J. Biol. 
5 Chem., 263: 9409-9417 (1988). The apparent Kd's are calculated using the method 
described by Apriletti (1995) and Apriletti (1988). The apparent Kd's are presented in 
TABLE 2. The apparent Kd's (App.Kd) are determined in the presence of the sample to 
be assayed, 1 nM [ l25 I]T 3 , and 50/xg/ml core histones, in buffer E (400 mM KC1, 200 
mM potassium phosphate, pH 8.0, 0.5 mM EDTA, 1 mM MgCl 2 , 10% glycerol, 1 mM 
10 DTT) in a volume of 0.21 ml. After incubation overnight at 4°C, 0.2 ml of the 

incubation mixture is loaded onto a Quick-Sep Sephadex G-25 column (2.7 x 0.9 cm, 1.7 
ml bed volume) equilibrated with buffer E. The excluded peak of protein-bound [ ,2S I)T 3 
is eluted with 1 ml of buffer E, collected in a test tube, and counted. Specific T, binding 
is calculated by subtracting nonspecific binding from total binding. 

15 
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20 Example 3 - Increased Nuclear Protein Coactivation by TR Ligands 

To test the ability of TR ligands to activate the binding of TR to the nuclear 
activation protein RIP-140 (a nuclear protein that can bind to nuclear receptors, such as 
the estrogen receptor), a TR ligand is liganded to TR and then incubated with RIP-140 as 
described in V. Cavailles, et ai, EMBO J., 14(15):3741- 3751 (1995), which is 

25 incorporated by reference herein. In this assay, 35 S -RIP-140 protein binds to liganded TR 
but not unliganded TR. Many TR 35 s ligands can activate RIP-140 binding as shown in 
TABLE 2. 
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Example 4 - TR ligand binding and TR activation in Cultured Cells 

To test TR activation of transcription in a cellular environment, TR ligands are 
assayed for their ability to activate a reporter gene, chloramphenicol transferase ("CAT"), 
which has a TR DNA binding sequence operatively linked to it. Either GC or L937 cells 
5 (available from the ATCC) can be used, respectively). In such assays, a TR ligand 
crosses the cell membrane, binds to the TR, and activates the TR, which in turn activates 
gene transcription of the CAT by binding the TR DNA binding region upstream of the 
CAT gene. The effective concentration for half maximal gene activation (EC^ is 
determined by assaying CAT gene activation at various concentrations as described herein 
10 and in the literature. The results of CAT gene activation experiments are shown in 
TABLE 2 



CAT GENE ACTIVATION ASSAYS 

Functional response to thyroid hormone O.S.S'-triiodo-L-thyronine, T 3 ) and TR 
15 ligands is assessed either in a rat pituitary cell line, GC cells, that contain endogenous 
thyroid hormone receptors (TRs) or U937 cells that contain exogenous TRs expressed as 
known in the art. GC cells are grown in 10-cm dishes in RPMI 1640 with 10% newborn 
bovine serum, 2 mM glutamine, 50 units/ml penicillin and 50 yig/ml streptomycin. For 
transfections, cells are trypsinized, resuspended in buffer (PBS, 0.1% glucose) and mixed 
20 with a TREtkCAT plasmid (10 mg) or phage in 0.5 mi buffer (15 ±5 million cells) and 
electroporated using a Bio-Rad gene pulser at 0.33 kvolts and 960 mF. The TREtkCAT 
plasmid contains two copies of a T 3 response element (AGGTCAcaggAGGTCA) cloned in 
the Hind III site of the pUC19 polylinker immediately upstream of a minimal (-32/ +45) 
thymidine kinase promoter linked to CAT (tkCAT) coding sequences. After 
25 electroporation, cells are pooled in growth medium (RPMI with 10% charcoal-treated, 
hormone stripped, newborn bovine serum), plated in 6-well dishes and treated with either 
ethanol or hormone. CAT activity is determined 24 hours later as described D. C. 
Leitman, R. C. J. Ribeiro, E. R. Mackow, J. D. Baxter, B. L. West, J. Biol. Chem. 
266, 9343 (1991), which is incorporated by reference herein. 
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Effect of TS-10 on the Transcriptional Regulation of the DR4-ALP 
Reporter Gene in the Presence or Absence of T3. 

Characteristics of the TRAF cells: TRAFal are CHO Kl cells stably transformed with 
an expression vector encoding the human thyroid hormone receptor a 1 and a DR4.ALP 
5 reporter vector; TRAFbl are CHO Kl cells stably transformed with an expression vector 
encoding the human thyroid hormone receptor 01 and a DR4-ALP reporter vector. 

Interpretation of the effect of compound TS-10 on the transcriptional regulation of 
the DR4-ALP reporter gene in the presence or absence of T3. 

10 

TRAFal reporter cells: TS-10 alone (open circles) induces a partial activation of the 
expression of the ALP reporter protein amounting to approximately 27% of the maximal 
effect by the natural thyroid hormone T3. In the presence of T3 (filled circles), TS-10 
has a weak antagonistic effect. The EC50 concentration for the agonistic effect of TS-10 
15 and the EC50 concentration for its T3 antagonistic effect, respectively, is indicated in 
FIG. 18. 

In FIG. 18, open and filled circles with dotted lines show the dose-dependent effect of 
TS-10/T3 on the toxicity marker (MTS/PMS), reduction of tetrazolium salt in the 
20 mitochondria, displayed on the right y-axis as optical density. There is no obvious toxic 
effect of TS-10 on the MTS-PMS marker but there is a clear effect on the morphology of 
the cells, as can be seen under the light microscope, at the highest concentration of TS-10 
(32 mM) both in the absence and presence of T3, respectively (not shown in the figure). 



25 TRAFbl reporter cells: TS-10 alone (open circles) induces a partial activation of the 
expression of the ALP reporter protein amounting to approximately 35% of the maximal 
effect by T3. The EC50 concentration for the agonistic effect of TS-10 is indicated in 
FIG. 19. In the presence of T3 (filled circles), TS-10 shows, if anything, a slight 
potentiation of the T3 effect on the expression of the ALP reporter protein. The T3 

30 inhibitory effect of TS-10 at its highest concentration used (32 mM) is a toxic effect 
rather than T3 antagonism. 
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In FIG. 19, open and filled circles with dotted lines show the dose-dependent effect of 
TS-10/T3 on the toxicity marker (MTS/PMS), reduction of tetrazolium salt in the 
mitochondria, displayed on the right y-axis as optical density. There is no obvious toxic 
effect of TS-10 on the MTS-PMS marker but a clear effect on the morphology of the cells 
5 can be observed, under the light microscope, at the highest concentration of TS-10 (32 
mM) both in the absence and presence of T3, respectively (not shown in the figure). 

HepG2 (HAF18) reporter cells: TS-10 alone (open circles) induces a partial activation 
of the expression of the ALP reporter protein amounting to slightly more than 50% of the 

10 maximal effect by T3. The EC50 concentration for the agonistic effect of TS-10 is 
indicated in FIG. 20. In the presence of T3 (filled circles), TS-10 shows no effect i.e. 
no T3 antagonism nor potentiation/additive effect to T3. Open and filled circles with 
dotted lines show the dose-dependent effect of TS-10/T3 on the toxicity marker 
(MTS/PMS), reduction of tetrazolium salt in the mitochondria, displayed on the right y- 

15 axis as optical density. There is no obvious toxic effect of TS-10 on the MTS/PMS 
marker or on the morphology of the cells, as can be observed using a light microscope, at 
any concentration of TS-10/T3 used. 

Example 5 - Comparisons of Human TR-a and Human TR-0 
20 Competition for [ m I}T 3 binding to TR LBD by T 3 and Triac 

The drug, triac, is a thyroid hormone agonist. Triac is 3,5,3'-triiodothyroacetic 

acid and is described in Jorgensen, Thyroid Hormones and Analogs in 6 Hormonal 

Proteins and Peptides, Thyroid Hormones at 150-151 (1978). Another compound that can 

be used in place of triac is 3,5-diiodo-3'-isopropylthyroacetic acid. Competition assays 
25 are performed to compare the displacement of [ 125 I]T 3 from binding with human TR-a 

LBD or human TR-0 LBD by unlabeled T 3 or triac. The results of such assays are 

depicted in FIG. 16. 

Standard binding reactions are prepared containing 1 nM [ ,2S I]T 3 , 30 fmol of 

human TR-a (empty symbols) or 0 (solid symbols), and various concentrations of 
30 competing unlabeled T 3 (circles) or triac (triangles). Assays are performed in duplicate. 
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Scatchard Analysis of [ 125 IjT 3 Binding to TR 

Human TR-a (left panel) or human TR-/3 (right panel) is assayed for T 3 binding in 

the presence of increasing concentrations of [ 125 1]T 3 . The apparent equilibrium 

dissociation constant (20 pM for a and 67 pM for fi) is calculated by linear regression 

5 analysis and is depicted in FIG. 17. 

3, 5-DlBROMO-4-(3%5'-DlISOPROPYL-4'-HYDROXYPHENOXY) BENZOIC ACID IS A 

TR-a Selective Synthetic Ligand. 



10 




15 3, 5-dibromo-4-(3\ 5 '-diisopropyl-4'-hydroxyphenoxy) benzoic acid (Compound 

11), the structure of which is drawn above, is assayed for binding to the two different 
isoforms of the TR, Tra and TR/?. Compound 11 exhibits an IC50 of 1.6 fiM for 
binding to TRa and an IC50 of 0.91 fiM for binding to TR0. Assays for determining 
selective binding to the TRa or TR/3 LBD can include reporter assays, as described 

20 herein. See also Hollenberg, et al., J. Biol. Chem., 270(24)14274-14280 (1995). 

Example 6 - Preparation and Purification of a TR-a LBD 

Rat TR-a LBD, residues Metl22 - Val410, is purified from E. coli ("LBD- 
122/410"). The expression vector encoding the rat TR-a LBD is freshly transfected into 

25 £. coli strain BL21(DE3) and grown at 22°C in a 50-liter fermenter using 2x LB 

medium. At an A^ of 2.5-3, IPTG is added to 0.5 mM and growth is continued for 3 h 
before harvesting. The bacterial pellet is quickly frozen in liquid nitrogen and stored at - 
70°C until processed. Extraction and purification steps are carried out at 4°C. The 
bacteria are thawed in extraction buffer (20MM Hepes, pH 8.-, 1 mM EDTA, 0.1% 

30 MTG, 0.1 mM PMSF, and 10% glycerol) at a ratio of 10 ml buffer/g bacteria. Bacteria 
are lysed by incubation for 15 min. with 0.2 mg/ml lysozyrhe and sonicated at maximum 
power while simultaneously homogenized with a Brinkmann homogenizer (Model PT 
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10/35 with generator PTA 35/2) until the solution loses its viscosity. After centrifugation 
for 10 min at 10,000 g. the supernatant is adjusted to 0.4 M KC1, treated with 0.6% PEI 
to precipitate fragmented DNA, and centrifuged for 10 min at 10,000 g. The rat TR-a 
LBD in the supernatant is then precipitated with 50% ammonium sulfate and centrifuged 
5 for 10 min at 10,000 g. The precipitate is resuspended with buffer B (20 mM Hepes, pH 
8.0, 1 mM EDTA, 1 mM DTT, 0.1 mM PMSF, 0.01% Lubrol, and 10% glycerol) to a 
final conductivity of 9 mS/cm (approx. 0.7 M ammonium sulfate) and centrifuged 1 h at 
100,000g. The supernatant is frozen in liquid nitrogen and stored at -70°C. 

The crude extract is thawed, bound with a tracer amount of [ ,25 I]T 3 , and loaded 
10 directly onto a phenyl-Toyopearl hydrophobic interaction column (2.6 x 18 cm, 95 ml bed 
volume) at 1.5 ml/min. The column is eluted with a 2-h gradient from 0.7 ammonium 
sulfate, no glycerol to no salt, 20% glycerol in buffer C (20 mM Hepes, pH 8.0, 0.5 mM 
EDTA, 1 mM DTT, 0.2 mM PMSF). The rat TR-a LBD prebound to tracer [ ,25 I)T 3 
(less than 0.005% of total rat TR-a LBD) is detected using a flow-through gamma 
15 emission detector, whereas unliganded rat TR-a LBD is assayed by postcolumn [ 125 I]T 3 
binding assays (described herein). 

The phenyl-Toyopearl unliganded rat TR-a LBD peak fractions are pooled, diluted 
with buffer B to a conductivity of 0.5 mS/cm (equivalent to approx. 20 mM ammonium 
sulfate), loaded onto a TSK-DEAE anion-exchange column (2 x 15 cm, 47 ml bed 
20 volume) at 4 ml/min, and eluted with a 60-min gradient from 50 to 200 mM NaCl in 
buffer B. 

The unliganded rat TR-a LBD peak fractions from TSK-DEAE are pooled, diluted 
twofold with buffer B, loaded at 0.75 ml/min on a TSK-heparin HPLC column (0.8 x 7.5 
cm, 3 ml bed volume), and eluted with a 50 to 400 mM NaCI gradient in buffer B. 

25 The pool of unliganded rat TR-a LBD peak fractions from the TSK-heparin 

column is adjusted to 0.7 M ammonium sulfate, loaded at 0.75 ml/min on a TSK-phenyl 
HPLC column (0.8 x 7.5 cm, 3 ml bed volume), and eluted with a 60-min gradient from 
0.7 M ammonium sulfate without glycerol to no salt with 20% glycerol in buffer C The 
fractions containing unliganded rat TR-a LBD are pooled and incubated with a five fold 

30 excess of hormone for 1 h, the salt concentration is adjusted to 0.7 M ammonium sulfate, 
and the sample is reloaded and chromatographed on the same column as described above. 
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Example 7 - Crystallization of Liganded TR-a LBD 

Material from a single LBD- 122/4 10 preparation is divided into batches, and 
quantitatively bound with one of the following ligands: Dimit, T 3 , or triac IpBr, 
(3,5dibromo-3'isopropylthyronine) for the final purification step. 

5 To maintain full saturation of rat TR-a LBD with a Hgand, and to prepare the 

complex for crystallization, the ligand-bound rat TR-a LBD is concentrated and desalted 
in an Amicon Centricon-10 microconcentrator (McGrath et al. Biotechniques, 7:246-247 
(1989), incorporated by reference herein), using 10 mM Hepes (pH 7.0), 3.0 mM DTT. 
and 1.0 nM to 10 nM ligand. 

10 Factorial crystallization screening trials (Jancarik & Kim, J. Appl. Crystallogr. 

24:409-411 (1991) incorporated by reference herein) are carried out for rat TR-a LBD 
bound to selected ligands using hanging-drop vapor diffusion at 17 °C (with 1 y\ protein 
solution, 1 Ail precipitant solution and a 0.5 ml reservoir using silanized coverslip: 
(McPherson, Preparation and Analysis of Protein Crystals (1982), incorporated by 

15 reference herein). Rat TR-a LBD is not stable at 4°C and is stored at -80°C, where it 
maintains its avidity for hormone and its crystallizability for approximately two to three 
months. These procedures are carried out as described in McGrath, M.E. et al., J. Mol. 
Biol. 237:236-239 (1994) (incorporated by reference).) Crystals are obtained in condition 
21 of the screening trials (Jancarik & Kim 1991) and conditions are then optimized. 

20 Wedge-shaped crystals are reproducibly obtained with hanging-drop vapor fusion at 22°C 
with 15% 2-methyl-2,4-pentanediol (MPD), 0.2 M ammonium acetate and 0.1 M sodium 
cacodylate (pH 6.7), 3 mM DTT, with 2 fi\ protein solution, 1 fil precipitant solution and 
a 0.6 ml reservoir using silanized coverslip, and with 8.7 mg/ml (Dimit), 5.5 mg/ml 
(IpBr 2 ), 5 mg/ml (triac), or 2.3 mg/ml (T 3 ) over a period of three days. Under these 

25 conditions, diffraction quality crystals (dimension 0.5 x 0.2 x 0.0075 mm 3 ) can be grown 
at ambient temperature (22 *C). The best crystals have a limiting dimension of 
approximately 100 urn and are obtained at a protein concentration between 2.3 and 8.7 
mg/ml in the presence of 3 mM DTT. The crystals are of the monoclinic space group 
C2, with one monomer in the asymmetric unit. 

30 
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Example 8 - Crystallization of Human TR-0 LBD Complexed with T 3 or 
Triac 

Human TR-/? LBD complexed with T 3 and human TR-0 LBD complexed with 
triac are purified according to the same procedures described above for the rat TR-a 
5 LBD, with the following modifications. 

The expression of human TR-/S LBD differs from the rat TR-a LBD in that the 
human TR-/3 LBD residues extend from the amino acid at position 716 through the amino 
acid at position 1022, according to the amino acid numbering scheme for the various 
nuclear receptor LBDs depicted in FIG. 3. FIG. 3 illustrates a numbering scheme 
10 applicable to all of the nuclear receptors listed as well as to any additional homologous 
nuclear receptors. The vertical lines on FIG. 3 at position 725 and at position 1025 
delineate the preferred minimum amino acid sequence necessary to obtain adequate 
binding of ligand. The amino acid sequence from position 716 to position 1022 according 
to the numbering scheme of FIG. 3 corresponds to the amino acid positions 202 to 461 
15 according to the conventional numbering of the amino acid sequence of human TR-0 
which is publicly available. Also, the human TR-/3 LBD is expressed with a histidine 
tag, as described in Crowe et al., Methods in Molecular Biology 31:371-387 (1994), 
incorporated by reference herein. 

The purification of human TR-0 LBD is the same as that described above for the 
20 rat TR-a LBD with the following exceptions. First, before the purification step using the 
hydrophobic interaction column, a step is added in which the expressed human TR-0 LBD 
is purified using a nickel NTA column (commercially available from Qiagen, Chatsworth, 
CA) according to manufacturer's instructions, and eluted with 200 mM imidazole. The 
second difference is that in the purification of the human TR-/3 LBD, the purification step 
25 using a heparin column is omitted. 

The crystallization of human TR-0 LBD bound to T 3 or triac is as follows. 
Crystals are obtained in condition 7 of the factorial screen using hanging drops as before 
at ambient temperature (22 °C) using the factorial crystallization screening trials of 
Jancarik & Kim (1991) and using the commercially available product from Hampton 
30 Research, Riverside). The following are optimum conditions: hexagonal bipyrimidal 
crystals are grown at 4°C for 2-3 days from hanging drops containing 1.0-1.2 M sodium 
acetate (pH unadjusted) and 0.1 M sodium cacodylate (pH 7.4), 3 mM DTT, with either 
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a 1 /xl protein solution, 1 ii\ precipitant solution or 2 jil protein solution, 1 /xl precipitant 
solution and a 0.6 ml reservoir using silanized coverslip, at a protein concentration of 7- 
10 mg/ml. The best crystals have a limiting dimension of 200 pm. 

The crystal system for human TR-/3 LBD bound to T 3 or triac is trigonal with the 
5 space group p3,21. The unit cell dimensions are cell length a = cell length b = 68.448 
angstroms, cell length c = 130.559 angstroms. The angles are a = 90°, (3 = 90°, 
gamma - 120°. 

Example 9 - Determination of Liganded TR-<* LBD Crystal Structure 

10 Data from each of three cocrystals (Rat TR-a LBD with Dimit, T3 and IpBr 2 ) is 

measured on a Mar area detector at Stanford Synchrotron Radiation Laboratory beamline 
7-1 (X = 1.08 angstroms) using 1.2* oscillations. 

Data from the T 3 cocrystal is measured with the b* axis approximately parallel 
with the spindle. The crystals are flash frozen at -178°C in a nitrogen gas stream with the 

15 MPD mother liquor serving as the cryosolvent. An orientation matrix for each crystal is 
determined using REFIX (Kabsch, W., /. Appl. Crystallogr. 26:795-800 (1993) 
incorporated by reference). Reflections are integrated with DENZO (commercially 
available from Molecular Structure Corp., The Woodlands, Texas), and are scaled with 
SCALEPACK (as described in Otwinowski, Z, Proceedings of the CCP4 Study Weekend: 

20 "Data Collection and Processing, " 56-62 (SERC Daresbury Laboratory, Warrington, UK 
1993) incorporated by reference). 

For the T 3 data set, Bijvoet pairs are kept separate, and are locally scaled using 
MADSYS (W. Hendrickson (Columbia University) and W. Weis (Stanford University)). 
Cocrystals prepared from the three isosteric ligands are isomorphous. MIR 

25 analysis is performed using programs from the CCP4 suite (Collaborative Computational 
Project, N.R. Acta Crystallogr. D50:760-763 (1994), incorporated by reference herein). 
Difference Pattersons is calculated for both T 3 and IpBr 3 , taking the Dimit cocrystal as the 
parent. The positions of the three iodine atoms in the T 3 difference Patterson are 
unambiguously determined from the Harker section of the density map as peaks of 11a 

30 above background. The positions for the two bromine atoms in the IpBr 2 cocrystals, are 
located independently, as peaks 8a above the noise level. Phases for the LBD-122/410 are 
calculated from the solution to the IpBr 2 difference Patterson, and are used to confirm the 
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location of the unique third iodine of the T 3 cocrystal. Halogen positions are refined with 
MLPHARE, including the anomalous contributions from the iodine atoms (Otwinowski, 
Z, Proceedings of the CCPR Study Weekend 80-86 (SERC Daresbury Laboratory, 
Warrington, UK 1991)). The MIRAS phases are improved through solvent 
5 flattening/histogram matching using DM (Cowtan, K., Joint CCP4 and ESF-EACBM 
Newsletter on Protein Crystallography 31: 34-38 (1994), incorporated by reference 
herein). 

A model of the LBD- 122/4 10 with Dimit bound is built with the program O from 
the solvent flattened MIRAS 2.5 angstrom electron density map (Jones et al., Acta 

10 Crystallogr. A 47:110-119 (1991), incorporated by reference herein). The initial model, 
without ligand, (Rcryst = 40.1%), is refined using least-squares protocols with XPLOR. 
The Dimit ligand is built into unambiguous Fo-Fc difference density during the following 
round. Subsequent refinement employs both least-squares and simulated annealing 
protocols with XPLOR (Brunger et al., Science 235:458-460 (1987), incorporated by 

15 reference herein). Individual atomic B-factors are refined isotropically. As defined in 
PROCHECK, all residues are in allowed main-chain torsion angle regions as described in 
Laskowski et al., J. Appl. Crystallogr. 26:283-291 (1993), incorporated by reference 
herein. The current model is missing 34 residues (Met I22 -Gln 156 ) at the N-terminus, and 5 
residues (Glu 406 -Val 410 ) at the C-terminus. 

20 In addition, the following residues are not modeled beyond Q3 due to poor 

density: 184, 186, 190, 198, 206, 209, 240, 301, 330, 337, 340, 343, 359, and 395. 
The average B-value for protein atoms is 34.5 A 2 . The final model consists of the LBD- 
122/410, residues Arg l57 -Ser m , Trp I85 -Gly 197 , Ser 199 -Asp 206 and Asp 208 -Phe 405 ; three 
cacodylate-modified cysteines: Cys 334 , Cys 380 and Cys 392 ; and 73 solvent molecules 

25 modeled as water (2003 atoms). 

*R, ym = 100 x £ hkl Eilli-H / Ehici Eil, 
tR* r = 100 x E hkI lF PH -F H l/ E hkl iFpl 

The occupancy for the two bromine sites is set to 35 electrons. The occupancies of the 
iodine sites are relative to this value. 
30 §Phasing power = (FH) / (e), where (FH) is the mean calculated heavy atom structure 
factor amplitude and (e) is the mean estimated lack of closure. 
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llRcullis = (e) / (iso), where <e> is the mean estimated lack of closure and (iso) is the 
isomorphous difference. 

iRcryst = 100 x £ hkl IF 0 -Fcl / £ hltl IF J where F 0 and F c are the observed and 
calculated structure factor amplitudes (for data F/a > 2). The Rfree was calculated using 
5 3% of the data, chosen randomly, and omitted from the refinement. 

§ Correlation coefficient = £ hkl flFj - DFJ) x (FJ - DFJ)/ Zw (PJ - ^J) 2 * Zm 

Example 10. Phasing of the rTRa LBD complex with Triac 

10 Due to the possible non-isomorphism of the rTRa LBD complex with Triac, a molecular 
replacement solution is determined using AM ORE (Navaza, J., Acta Crystallographica 
Section A-Fundamentals of Crystallography 50:157-63 (1994) from a starting model 
consisting of rTRa LBD complex with T 3 , but with the ligand, all water molecules, and 
the following residues omitted: Asn 179, Arg228, Arg262, Arg266, and Ser 277. Strong 

15 peaks are obtained in both the rotation and translation searches, with no significant (> 
0.5 times the top peak) false solutions observed (Table 3). Strong positive density present 
in both the anomalous and conventional difference Fourier maps confirm the solution. 
Maps are calculated using sigma-A weighted coefficients output by REFMAC 
(Murshudov, et al. "Application of Maximum Likelihood Refinements," in the 

20 Refinement of Protein Structures, Proceedings of Daresbury Study Weekend (1996)) after 
15 cycles of maximum likelihood refinement. Triac, the omitted residues, and water 
molecules 503, 504, 534 (following the numbering convention for the TR complex with 
T3) are built into the resulting difference density using O (Jones et. al.); the 
conformations of these residues are further confirmed in a simulated-annealing omit map 

25 (Brunger et. al.). The complete model is then refined using positional least-squares, 
simulated annealing, and restrained, grouped B factor refinement in XPLOR to an Rcryst 
of 23.6% and an Rfree of 24.1% 

Example 11. Connecting QSAR with Structure in the Thyroid Hormone 
30 Receptor 

The conclusions of classic thyroid hormone receptor quantitative structure-activity 
relationships may be summarized as follows: 
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1) the R/-hydroxyl group functions as a hydrogen bond donor; 

2) the amino-propionic acid interacts electrostatically through the carboxylate 
anion with a positively charged residue from the receptor; 

3) the preferences of R 3 /R 5 substituent are I > Br > Me > > H; 

5 4) the preferences of the R 3 ' -substituent are Ipr > I > Br > Me > > H. 

The structure of the thyroid hormone receptor ligand binding domain complexed with the 
agonists 3,5,3'-triiodothyronine (T 3 ), 3,5-dibromo-3'-isopropylthyronine (IpBr 2 ), 3,5- 
dimethyI-3'-isopropylthyronine (Dimit), and 3,5,3*-triiodothyroacetic acid (Triac), as 
provided herein, permits: 
10 1) the identification of receptor determinants of binding at the level of the 

hydrogen bond; 

2) the association of these determinants with the predictions of classic thyroid 
hormone receptor QSAR; and 

3) prediction as to which determinants of binding are rigid, and which are 
15 flexible, for both the ligand and the receptor. 

This classification for the agonists of the type (Rj= amino-propionic, acetic acid; 
R 3 ,R 5 = I,Br,Me; R 3 '=Ipr,I) is given below (for the representative ligand T 3 ); 



F = Fiducial (always satisfied) 
20 A = Adjustable 



25 




i i 

A A 



30 
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Based upon the methods and data described herein, the following is an 
embodiment of the computational methods of the invention, which permit design of 
nuclear receptor ligands based upon interactions between the structure of the amino acid 
residues of the receptor LBD and the four different ligands described herein. The small 
5 molecule structures for the ligands can be obtained from Cambridge Structural Database 
(CSD), and three dimensional models can be constructed using the methods described 
throughout the specification. The following are factors to consider in designing synthetic 
ligands: 

1) Histidine 381 acts as a hydrogen bond acceptor for the R,' hydroxyl, with the 

10 optimal tautomer maintained by water molecules. See FIG. 23 and FIG 24. Histidine is 
the only hydrophilic residue in this hydrophobic pocket that surrounds the IV substituent. 
Histidine can be either a hydrogen bond acceptor or donor, depending on its tautomeric 
state. It is preferably a hydrogen bond donor , but can tolerate being a hydrogen bond 
acceptor, as for example, when there is a methoxy at the R/ position of the ligand; 

15 2) Arginines 228, 262, and 266 interact directly and through water-mediated 

hydrogen bonds with the R r substituent, with the electrostatic interaction provided by 
Arginine 266 (as in the Triac complex). This polar pocket is illustrated by FIG. 23 - 
FIG. 25. FIG- 23 depicts T 3 in the TRa ligand binding cavity, where T3's amino- 
propionic Rl- substituent interacts with Arg 228, H0H502, H9H503 and H0H5O4 via 

20 hydrogen bonds. FIG. 24 depicts triac in the ligand binding cavity, with its -C00H R x 
substituent in the polar pocket. In FIG. 24, Arg 228 no longer shares a hydrogen bond 
with the ligand, but the -COOH R { substituent forms hydrogen bonds with Arg 266. 
FIG. 25 superimposes T 3 and triac in the ligand binding cavity and shows several 
positionally unchanged amino acids and water molecules, and selected changed interacting 

25 amino acids and water molecules. The three figures illustrate pans of the polar pocket 
that can change and those parts that do not move upon binding of different ligands. For 
example, the Arg 262 at the top of the polar pocket does not move, even when the R, 
substituent has changed from a -COOH to an aminopropionic acid group. However, the 
other two Arginines, Arg 228 and Arg 266, demonstrate flexibility in the polar pocket to 

30 respond to the change in the size or chemical naure of the R t substituent.; 

3) Inner and outer pockets for the R 3 /R 5 substituents are formed by Ser260, 
Ala263, Ile299; and Phe 218, Ile221, Ile222, respectively. See FIGS. 21 and 22. The 
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inner pocket is filled by either the R 3 or the R 5 substituent, regardless of the size of the 
substituent, and may act as a binding determinant by positioning the ligand in the 
receptor. Optimally, the inner pocket amino acids interact with an R3 or R5 substituent 
that is no larger than an iodo group. If the inner pocket is filled by the R 3 substituent, 
5 then the outer pocket interacts with the R 5 substituent and vice versa. The outer pocket 
can adjust to the size of its substituent through main chain motion centered at the break in 
helix 3 (Lys220-Iie221), suggesting that the bending of H3, and motion of the N-terminal 
portion of H3, may represent a conformational change induced on ligand binding. The 
outer pocket has greater flexibility than does the inner pocket in terms of accommodating 

10 a larger substituent group. 

4) A pocket for the R 3 *-substituent is formed by Phe 215, Gly290, Met388. 
The pocket is incompletely filled by the R 3 *-iodo substituent, and accommodates the 
slightly larger 3'-isopropyl substituent by movement of the flexible Met388 side chain and 
the H7/H8 loop. This pocket can accommodate R 3 * substituents that are even larger than 

15 isopropyl, for example, a phenyl group. 

The above information will facilitate the design of high affinity agonists and 
antagonists by improving automated QSAR methodologies and informing manual 
modeling of pharmaceutical lead compounds. For example, the inclusion of discrete 
water molecules provides a complete description of hydrogen bonding in the polar pocket 

20 for use with pharmacophore development: also, the identification of mobile and immobile 
residues within the receptor suggests physically reasonable constraints for use in 
molecular mechanics/dynamics calculations. 

Example 12. Design of an Increased Affinity Ligand 
25 Direct interaction between the receptor and the ligand is limited in the polar 

pocket, which interacts with the R, substituent. While the lack of complementarity may 
contain implications for biological regulation, it also provides an opportunity for 
increasing affinity by optimizing the interaction between the amino acids of the polar 
pocket and the R, substituent of a synthetic ligand. The structure of the receptor-ligand 
30 interactions described herein enables design of an increased affinity synthetic ligand 
having two complementary modifications: 
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1) Remove the positively charged amine. The strongly positive electrostatic 
potential predicted for the polar pocket suggests that the positively charged amine of the 
aminopropionic acid R, substituent may be detrimental to binding. Suitable groups for 
substitution are suggested by the nature of nearby hydrogen bond partners: for example, 

5 Thr 275 O or Ser 277 N. See e.g. Tables in Appendix 2. For example, any any 

negatively charged substituent would be compatible for interacting with the amino acids of 
the polar pocket, including carboxylates, carbonyl, phosphonates, and sulfates, comprising 
0 to 4 carbons. Another example of an R, substitution is an oxamic acid that replaces the 
amine of the naturally occurring ligand with one or more carbonyl groups. 

10 2) Incorporate hydrogen bond acceptor and donor groups into the R,- 

substituent to provide broader interactions with the polar pocket scaffold. Such hydrogen 
bond acceptor and donor groups incorporated into the Rl -substituent will allow 
interactions that would otherwise occur with water molecules in the polar pocket. 
Specific waters include HOH 504 (hydrogen bonds with Ala 225 O and Arg 262 NH); 

15 and HOH 503 hydrogen bonds with Asn 179 OD1, Ala 180 N), both of which are present 
in all four complexes (TR LBD complexed with T3, TR LBC complexed with IpBr 2 , TR 
LBD complexed with Dimit and TR LBD complexed with Triac). Analysis of the 
hydrogen bonding network in the polar pocket suggests replacement of HOH 504 with a 
hydrogen bond acceptor, and HOH 503 with an hydrogen bond donor (although the 

20 chemical nature of asparagine probably permits flexibility at this site). Thus, 

incorporating a hydrogen bond acceptor in an Rl substituent that could take the place of 
the HOH504 or incorporating a hydrogen bond acceptor in an Rl substituent that could 
positionally replace the HOH503, or a combination thereof, are methods of designing 
novel synthetic TR ligands. 

25 These two design approaches can be used separately or in combination to 

design synthetic ligands, including those in Table 4 (below). 

A corollary to this approach is to design specific interactions to the residues 
Arg262 and Asn 179. The goal is to build in interactions to these residues by designing 
ligands that have R, substituents that form hydrogen bonds with water molecules or 

30 charged residues in the polar pocket. 



WO 97/21993 



68. 

Table 4: Synthetic TR Ligands 
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s"R R* 6 Rs Re 

3'R R*2 R3 RjJ 



R1 R2 R3 R5 R6 X R'2 R'3 R'4 R'5 R'6 



10 C02H 



Me Me H 0 H Me OH Me 



CH2C02H II S - El SH Et 

CH2CH2C02H Br Br nPr NH2 nPr 

CH2CH(NH2)C02H CI CI iPr iPr 

OCH2C02H Et Et Ph nB u 

15 OCH2CH2C02H OH OH I n Pen 

NHCH2C02H NH2 NH2 Br nHex 

NHCH2CH2C02H SH SH CI Ph 



CH2COCOC02H 



hetero 
cycle 



20 NHCOCOC02H ary , 
COC02H 
CF2C02H 
COCH2C02H 



25 Any combination of the above substituents in the biphenyl ether scaffold structure shown above may result in a 
potentially pharmacologically useful ligand for the thyroid hormone receptor. These novel ligands may be 
antagonists of the thyroid receptor. 
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69. 

cgy" for designing synthetic ligands using the computational methods described 
is summed below: 



A = Hydrogen Bond Acceptor 

D = Hydrogen Bond Donor 

O = -OH, -CO 

RIO can be -OH, -CO 

R20 can be -CO 

R30 can be -COOH, -CONH2 

See also Table of synthetic TR Ligands 




For example. 
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10 



15 



20 



25 



30 



Dimil 

Data collection 

Cell dimensions 

a (A) 117.16 

b (A) 80.52 

c (A) 63.21 

P C) 120.58 

Resolution (A) 2.2 

Obs. Reflections. 57031 
(no.) 

Unique Reflections. 22327 
(no.) 

Completeness, (%) 87.0 

*R» r n (%) 3.9 

Phasing (15.0 - 2.5A) 

No. of sites 
^Occupancy 
(Anomalous) 

§F H /E 

centric (acentric) 
15.0-5.0 A 
5.0-3.0 A 
3.0-2.5 A 

15.0-5.0 A 

5.0-3.0 A 

3.0-2.5 A 

Mean figure of merit 0.62 

MR Phasing 
(10-3.5A) 



35 



Rotation Search: 
Eyler Angles (*) 



T3 



117.19 

80.20 

63.23 

120.60 

2.0 

64424 

21023 

82.4 
3.5 

19.6 
3 

44.6 (19.8) 
50.2 (23.7) 
39.2 (22.3) 



3.67 (4.61) 
2.23 (2.75) 
1.64 (1.99) 

33 
45 
60 



IpBr ; 



117.18 

80.12 

63.13 

120.69 

2.1 

66877 

23966 

93.7 
4.5 

11.6 
2 

35.0 
35.0 



2.25 (3.09) 
1.25 (1.85) 
1.15 (1.57) 

44 

63 
65 



Triac 



118.19 

81.37 

63.73 

121.00 

2.45 

83573 

18453 

96.0 
7.5 



G, = 309.37 
= 48.96 
0 3 = 127.28 



40 



§ correlation 
coefficient 

Translation Search. 

Fractional 
coordinates 



34.3 

x = 0.1571 
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y - 0.000 



z - 0.3421 



§ correlation 
Coefficient 

5 'R factor 



65.8 



31.2 



Refinement 
Resolution (A) 



15.0-2.2 



5.0-2.0 



15.0 - 2.2 



25-2.5 




22.1 



21.4 



23.6 



22.7 



24.0 



22.4 



24.1 



10 



All publications and patent applications mentioned in this specification are herein 
incorporated by reference to the same extent as if each individual publication or patent 
application was specifically and individually indicated to be incorporated by reference. 
15 The nuclear receptor ligands, particularly the TR ligands, of these references are herein 
incorporated by reference and can be optionally excluded from the claimed compounds 
with a proviso. 

Headings and subheadings are presented only for the convenience of the reader 
and should not be used to construe the meaning of terms used within such headings and 
20 subheadings. 

The invention now being fully described, it will be apparent to one of ordinary 
skill in the art that many changes and modifications can be made thereto without 
departing from the spirit or scope of the appended claims. 
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.c-rawcns or Oimn v.r,n rhyro.a Homwr.i .neceptor Afr.ino Acids 



Dimit 
Atom 



Amino Acid 
In full length a 



Amino Acid 
Atom 



Distance 



C16 

C16 

CI9 

C18 

C18 

CI9 

C18 

C16 

CI9 

ce- 
cs 

..C10 



C12 
C13 

03 
C10 
C10 

cn 



215-PHE 
215-PHE 
218-PHE 
218-PHE 
218-PHE 
218-PHE 

218- PHE 

219- THR 
22MLE 
222-ILE •:. : 
222-ILE 
222-ILE: 
222-iLE 

■ -• 222-ILE t- ■ 
225-ALA 
225iALA ; 
228-ARG 
- . : ; 228-ARG;/;;v.: 
228-ARG 
:;.:228 : ARGS4 
258-MET 
-vv^2^ETiS 

258- MET 

■^/ : :-258-MEM,> 

259- MET 



^::SiWvP«i ::: €! 1- ' :25^ME« 



C15 
C1 
C13 

£:1E:C^;::;:j 

C17 
04 



259-MET 

259^METE:-i:e 
2S9-MET 

;259-MET|i:S 
259-MET 
:262-ARG: : i^^ 
262-ARG 

:«^argM: 

262-ARG 



3.98 
3.88 
3.69 
3.89 
3.92 
4.13 
3.77 
3.68 
4.11 
4.18 
3.72 
3.533.- 
3.85' 
' 4.13%j 
3.64 

4:02:iS,. :f £ " 
3.96 

3.36fe^/:S: 

3.58 

;; :3^91PS:JS1 
3.88 



3.42 

' 3.7' 
4^0 

4.09 



: ^smm i»» r c0 --• - 

^:^.3 ; 7OT- fefe 



C17 



03 



C7 
C9 

£ ^£4 ;cn>: 

N1 
C13 

C1S 



^•-•v ; \:4.03^E'3 : ": 

3.58. ........ 

3.85 . 
3.61 



263-ALA 



N1 276-LEU 



276-LEU fer V 
278-LEU 
278-LEU 

:IF':J':276-LEUSi 

276- LEU 

2n : sERvxs 

277- SER 



CA 
NH1 

cf 

C02 
CD2 

C02 

N 



3.69 

3.93 

-;;-3.82K 

" 3.51 

-"•".3. 
4.05 

4.09 
3.95 

:"4.'14-v^ ; 
3.79 
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10 



15 



20 



25 



30 



35 



40 



45 



Dimit 


Ammo Acid 






Atom 


in full length a 


Atom 


Distance A 


C17 


277-SER 


N 


3.69 


N1 


277-SER 


N 


3.30 


03 


277-SER 


N 


3.19 


C17 


277-SER 


CA 


3.92 


03 


277-SER 


CA 


3.35 


C13 


277-SER 


OG 


3.92 


C7 


287-LEU 


C02 


3.90 


C18 


290-GLY 


C 


4.04 


C18 


290-GLY 


0 


3.54 


C18 


291-GLY 


CA 


4.04 


C18 


292-LEU 


N 


4.20 


C2 


292-LEU 


CG 


4.18 


C4 


292-LEU 


CG 


3.88 


C6 


292-LEU 


CG 


4.01 


C2 


292-LEU 


CD1 


3.88 


C4 


292-LEU 


CD1 


4.02 


02 


292-LEU 


CD1 


4.07 


C4 


292-LEU 


CD2 


4.05 


C6 


292-LEU 


CD2 


3.72 


C8 


292-LEU 


CD2 


3.69 


CIO 


292-LEU 


CD2 


3.98 


01 


292-LEU 


CD2 


4.16 


C20 


299-ILE 


CD1 


3.87 


C8 


381-HIS 


CD2 


3.90 


C10 


381-HIS 


CD2 


3.84 


01 


381-HIS * 


G02 


3.40 


01 - 


381-HIS 


CE1 


3.72 


C8 


381-HIS 


NE2 


3.47 


C10 


381-HIS 


NE2 


3.51 


01 


381-HIS 


NE2 


2.64 


C8 


388-MET 


CE 


3.90 ! 


C8 


401-PHE 


CE1 


4.19 


01 


401-PHE 


CE1 


3,37 


C16 


401 -PHE 


CZ 


3.97 


01 


401-PHE 


CZ 


3.28 


C17 


3-H 2 0 


01 


3.77 


03 


502 


01 


3.13 


04 


3-H 2 0 


U 1 




C15 


2-H 2 0 


01 


4.04 


C17 


2-H 2 0 


01 


3.84 


N1 


2-H 2 0 


01 


3.35 


03 


503 


01 


2.56 


C17 


4-H 2 0 


01 


3.92 


04 


504 


01 


2.72 



Legend to Table XX. The table lists the interactions with Dimit (DMT). The column headings 
50 are as follows: 

#1 The atom of Dimit that interacts with the amino acid of the receptor. These are also 
numbered in figure X. 

U2 The amino acid in the full length rTRa that interacts with the ligand. 
#3 The name of the atom in the amino acid (standard nomenclature) where the interaction 
55 occurs. 

#4 The distance in A between Dimit and the protein atom. 
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i r lac 


Amino Acid 


Amino Acid 




Atom 


in full length q 


Atom 


Distance A 


1 1 


2 1 8-PHE 


0 


3.52 


II 


221-ILE 


CD1 


4.16 


■ i 
n 


222-ILE 


CA 


4.15 


1 1 


222-ILE 


CB 


4.03 


n 


221-ILE 


CG1 


3.92 


Lo 


222-ILE 


CD1 


4.12 


C 1 0 


222-ILE 


CD1 


3.77 


Liz 


222-ILE 


CD1 


3.79 


pii 
11 J 


225-ALA 


CB 


4.17 




225-ALA 


CB 


3.86 


C 1 0 


258-MET 


SD 


3.45 


C12 


258-MET 


SD 


3.73 


C10 


258-MET 


CE 


3.66 


C12 


256-MET 


CE 


3.77 


Cll 


256-MET 


C 


3.68 


Ctl 


256-MET 


0 


3.24 


CI 


259-MET 


0 


3.93 


i C1 1 


259-MET 


0 


3.24 


03 


259-MET 


0 


4.09 


C1 


259-MET 


CB 


3.89 


C13 


259-MET 


0 


3.74 


C14 


259-MET 


0 


3.96 


C1 


259-MET 


CB 


3.89 


Cll 


259-MET 


CB 


3.68 


C13 


259-MET 


CB 


4.01 


C14 


262-ARG 


CB 


4.07 


04 


262-ARG 


CB 


3.60 


C17 


282-ARG 


CD 




04 


262-ARG i 


CD 




03 


263-ALA 


N 


3.79 


03 


263-ALA 


CA 


3.67 


03 


283-ALA 


CB 


3.49 


C11 


107-ALA 


CB 


4.00 


03 


107-ALA 


CB 


3.49 


03 


266-ARG 


NH1 


3.00 


03 


275-THR 


0 


3.20 


03 


278-LEU 


CA 


3.11 


N1 


276-LEU 


C 


3.52 


03 


1 20-LEU 


N 


4.04 


C14 


1 20-LEU 


CA 


3.98 


03 


1 20-LEU 


CA 


3.11 


C14 


1 20-LEU 


c 




03 


1 20-LEU 


CB 


3.95 


02 


276-LEU 


CD1 


4.03 


11 


276-LEU 


CD1 


4.10 


C7 


276-LEU 


CD2 


3.84 


C9 


276-LEU 


CD2 


3.73 


Cll 


276-LEU 


CD2 


4.06 


02 


276-LEU 


CD2 


4.10 
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15 



20 



25 



30 



Tnac 
Atom 


Amino Acid 
in full length a 


Ammo Acid 
Atom 


======= 

uisiance m 


C13 


277-SER 


N 




C14 


277-SER 


N 


o, 1 o 


04 


277-SER 


N 


o./o 


03 


277-SER 


N 


o.Ub 


C14 


277-SER 


CA 


/O 


04 


277-SER 


CA 




C3 


277rSER 


OG 


1 R7 


12 


290-GLY 


0 


*3 c "J 
O.D / 


12 


292-LEU 


CG 




C4 


292-LEU 


CG 


3.9b 


C6 


292-LEU 


CG 


3.6b 


C2 


292-LEU 


CDV 


4.1 I 


C4 


292-LEU 


CD1 


3.85 


12 


292-LEU 


CD2 


3.98 


C4 


292-LEU 


CD2 


4.1 1 


C6 


292-LEU 


CD2 


*> A A 

3.44 


C8 


292-LEU 


CD2 


3.2o 


CIO 


292-LEU 


CD2 


3.88 


01 


292-LEU 


CD2 


3.35 


13 


299-ILE 


CD1 


3.77 


C8 


381-HIS 


CD2 


3.87 


C10 


381-HIS 


CD2 


3.90 


01 


381-HIS 


G02 


3.20 


01 


381-HIS 


CE1 


3.81 


C8 


381-HIS 


NE2 


3.51 


C10 


381-HIS 


NE2 


3.52 


01 


381-HIS 


NE2 


2.64 


01 


388-MET 


CE 


4.03 


01 


401 -PHE 


CE1 


3.86 


01 


401 -PHE 


CZ 


3.70 



Legend to Table XX. The table lists the interactions with triac. The column headings are as 
follows: 

#1 The atom of triac that interacts with the amino acid of the receptor. These are also 
35 numbered in figure X. 

02 The amino acid in the full length rTRa that interacts with the ligand. 

#3 The name of the atom in the amino acid (standard nomenclature) where the interaction 



occurs. 

#4 The distance in A between triac and the protein atom. 
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IpBR; Atom 


| Amino Acid 
j in full length a 


Amino Acid 
Atom 


Distance 
A 


C16 


215-PHE 


CD1 


4.01 


C16 


215-PHE 


CE1 


3.78 


BR1 


218-PHE 


0 


3.24 


BR1 


218-PHE 


c 


3.98 


C16 


218-PHE 


CB 


3.81 


C18 


218-PHE 


CB 


3 92 


BR1 


218-PHE 


CB 


4.08 


C18 


218-PHE 


CD2 




C18 


219-THR 


CG2 




BR1 


221-ILE 


CG1 


o.o \ 


BR1 


221-ILE 


CD1 


A, H7 


C6 


222-ILE 


CD1 




C8 


222-ILE 


CD1 




C10 


222-ILE 




O.DU 


C12 


222-ILE 


CD1 




01 


222-ILE 


CD1 


A OR 


C13 


225-ALA 


CB 


3 76 


04 


225-ALA 


CB 


*♦ .W 1 


04 


228-ARG 






C17 


228-ARG 


NH2 




03 


228-ARG 


NH2 


3 43 


04 


228-ARG 


NH2 


2.79 


C10 


256-MET 


SD 


3.65 


C12 


256-MET 


SD 


3.71 


C10 


256-MET 


CE 


3 90 


C12 


256-MET 


CE 


3.75 


BR2 


256-MET 


CE 


4.03 


C11 


259-MET 


c 


3.98 


C11 


259-MET 


0 


3.52 


C15 


259-MET 


0 


3.44 


N1 


259-MET 


0 


3.76 


C11 


259-MET 


CB 


3.87 


C15 


262-ARG 


CB 


4.03 
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ipBR 2 Atom 


Amino Acid 
in full length a 


Amino Acid 
Atom 


Oistance 

A 


C17 


262-ARG 


CB 


3.56 


03 


262-ARG 


CB 


3.55 


04 


262-ARG 


CB 


3.91 


C17 


262-ARG 


CD 


4.09 


04 


262-ARG 


CD 


3.71 


N1 


263-ALA 


N 


3.61 


N1 


263-ALA 


CA 


3.59 


N1 


263-ALA 


CB 


3.54 


03 


268-ARG 


NH 1 


3.93 


N1 


275-THR 


0 


3.43 


N1 


276-LEU 


CA 


3.46 


N1 


276-LEU 


C 


3.83 


C5 


276-LEU 


CD! 


4.02 


C7 


276-LEU 


CD2 


4.00 


C9 


276-LEU 


CD2 


3.81 


C11 


276-LEU 


CD 2 


3.91 


C13 


277-SER 


N 


3.79 


C15 


277-SER 


N 


3.63 


C17 


277-SER 


N 


3.70 


N1 


277-SER 


N 


3.17 


03 


277-SER 


N 


3.37 


C17 


277-SER 


CA 


3.89 


03 


277-SER 


CA 


3.43 


C13 


277-SER 


OG 


3.66 


02 


287-LEU 


CD1 


4.05 


C18 


290-GLY 


C 


4.04 


C18 


290-GLY 


0 


3.48 


C18 


291 -GLY 


CA 


4.02 


C4 


292-LEU 


CG 


3.89 


C6 


292-LEU 


CG 


a no 


C2 


292-LEU 


CD1 


3.79 


C4 


292-LEU 


CD1 


3.96 


02 


292-LEU 


CD1 


3.97 
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10 



15 



20 



25 



l n DD Ai*m 


Amino Acic 


Amino Acid 




Distance 




in full length o 


Atom 


A 


C4 


292-LEU 


CD2 


4.07 


C6 


292-LEU 


CD2 


3.75 


C8 


292-LEU 


CD2 


3.67 


CIO 


292-LEU 


CD2 


3.92 


BR2 


299-ILE 


CD1 


3.68 


C8 


361-HIS 


CD2 


3.67 


C10 


381-HIS 


CD2 


3.92 


01 


381-HIS 


G02 


3.50 


01 


381-HIS 


CE1 


3.62 


C8 


381-HIS 


NE2 


3.36 


CIO 


381-HIS 


NE2 


3.34 


01 


381-HIS 


NE2 


2.62 


C8 


401 -PHE 


CE1 


4.02 


01 


401 -PHE 


CE1 


3.19 


C16 


401 -PHE | 


CZ 


4.03 


01 


401 -PHE 


CZ 


3.06 


C17 


3-H 2 0 


01 




03 


502H 2 O 


01 


3.40 


04 


3-H 2 0 


01 




C15 


2-HjO 


01 




C17 


2-H 2 0 


01 




N1 


502H 2 0 


01 


3.12 


03 


503H 2 O 


01 


2.27 


C17 


4-H 2 0 


01 




04 


504H 2 0 


01 


1.78 



Legend to Table XX. The table lists the interactions with lpBr2. The column headings are as 
follows: 

#1 The atom of JpBr2 that interacts with the amino acid of the receptor. These are also 
30 numbered in figure X. 

#2 The amino acid in the full length rTRo that interacts with the iigand. 

#3 The name of the atom in the amino acid (standard nomenclature} where the interaction 

occurs. 

#4 The distance in A between lpBr2 and the protein atom. 
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T3 Atom 


Amino Acid 


Amino Acid 


Distance 




in full length o 


Atom 


A 


12 


215-PHE 


CD1 


4.08 


11 


218-PHE 


0 


3.19 


11 


218-PHE 


CB 


3.99 


C4 


218-PHE 


CB 


4.04 


11 


218-PHE 


CB 


3.99 


n 


221-ILE 


CG1 


4.01 


C8 


222-ILE 


CD1 


3.99 


CIO 


222-ILE 


CD1 


3.99 


C12 


222-ILE 


CD1 


3.57 


01 


222-ILE 


CD1 


3.68 


C13 


225-ALA 


CB 


3.66 


C3 


225-ALA 


CB 


4.04 


04 


228-ARG 


NH1 


3.23 


04 


228-ARG 


CZ 


3.45 


C15 


228-ARG 


NH2 


3.54 


03 


228-ARG 


NH2 


3.90 


04 


228-ARG 


NH2 


2.86 


C10 


256-MET 


SD 


3.73 


C12 


256-MET 


SO 


3.90 


C10 


256-MET 


CE 


3.97 


C12 


256-MET 


CE 


3.92 


C11 


259-MET 


C 


3.95 


C11 


259-MET 


0 


3.59 


C14 


259-MET 


0 


3.51 


N1 


259-MET 


0 


3.88 


CI 


259-MET 


CB 


4.06 


C11 


259-MET 


CB 


3.77 


C13 


259-MET 


CB 


3.96 


C15 


262-ARG 


CB 


3.61 


C14 


262-ARG 


CB 


4.02 


03 


262-ARG 


CB 


3.65 


04 


262-ARG 


CB 


3.92 


04 


282-ARG 


CD 


3.72 



APPENDIX 2 



WO 97/21993 



84. 



PCT/US96/20778 



15 



20 



25 



30 





1 O MiOiTl 


Amino Acid 


Amino Acid 




Distance 






in full length a 


Atom 


A, 




Nl 


263-ALA 


N 


J.O 1 




N1 


263-ALA 


CA 


o.o 1 




Nl 


263-ALA 


CB 


J.Oo 


5 


Nl 


275-THR 


o 




Nl 


276-LEU 


CA 


T TO 

O.JO 




Nl 


276-LEU 


p 


3.73 




C5 


276-LEU 


pm 


4.00 




C7 


276-LEU 


pm 


4.05 




C7 


276-LEU 


CD2 


3.80 


) 


C9 


276-LEU 




3.70 




C11 


276-LEU 




4.01 




C14 


277-SER 


M 

N 


3.67 




C15 


277-SER 


M 


3.62 




04 


228-ARG 


k|U1 


3.23 




Nl 


277-SER 


IN 


3.07 




03 


277-SSER 


M 
IN 


3.24 




C15 


277-SER 


PA 


3.77 




03 


277-SER 


PA 


3.34 




C13 


277-SER 


wo 


3.92 




!2 


290-GLY 


o 


3.50 




C4 


292-LEU 


pn 


O ft c 

3.95 




. C8 


292-LEU 


pr; 


3.83 




C2 


292-LEU 




4.07 




C4 


292-LEU 


pm 


3.99 




C4 


292-LEU 


PD9 


4.09 




C6 


292-LEU 




3.58 




C8 


292-LEU 


CD2 


J. 50 




CIO 


292-LEU 


CD2 


3.96 




01 


292-LEU 


CD2 


3.71 




13 


299-ILE 


CD1 


3.74 




C8 


381 -HIS 


CD2 


3.94 




C10. 


381 -HIS 


CD2 


3.97 




01 


381-HIS 


CD2 


3.39 
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T3 Atom 


Amino Acid 


Amino Acid 


Distance 




in full length a 


Atom 


A 


01 


381-HIS 


CD1 


3.82 


C8 


381-HIS 


NE2 


3.47 


C10 


381-HIS 


NE2 


3.55 


01 


381-HIS 


NE2 


2.70 


01 


388-MET 


CE 


3.88 


01 


401 -PHE 


CE1 


3.52 


01 


401 -PHE 


CZ 


3.32 


03 


502 


01 


2.51 


04 


3-H 2 0 


01 




N1 


2-H 2 0 


01 




03 


503 


01 


2.81 


04 


504 


01 


2.73 



Legend to Table XX. The table lists the interactions with T3. The column headings are as 
15 follows:' 

#1 The atom of T3 that interacts with the amino acid of the receptor. These are also 
numbered in figure X. 

#2 The amino acid in the full length rTRa that interacts with the ligand. 
#3 The name of the atom in the amino acid (standard nomenclature) where the interaction 
20 occurs. 

#4 The distance in A between T3 and the protein atom. 
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Coordination Structure of Thyroid Hormone Receptor and Triac 



Coordination Structure 
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Coordination Structure of Thyroid Hormone Receptor and IpB:2 
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Coordination Structure of Thyroid Hormone Receptor and T3 
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90. 



WHAT IS CLAIMED IS: 



5 



A crystal of a TR LBD comprising: 

1) a TR LBD and a TR LBD ligand or 

2) a TR LBD without a TR LBD ligand; 

wherein said crystal diffracts with at least a 2.0 to 3.0 angstrom resolution and has 



a crystal stability within 5% of its unit cell dimensions. 

2. The crystal of claim 1 , wherein said TR LBD has at least 200 amino acids. 

10 3. The crystal of claim 2, wherein said TR LBD is from a TR protein selected from 
the group consisting of TR amino acid sequence 122 to 410 of rat TR-a, 157 to 410 rat 
TR-a and 202 to 461 of human TR-/?. 

4. The crystal of claim 2, wherein said TR LBD ligand is a compound being of a 
15 formula: 



20 




25 



wherein R, is 

-0-CH 2 CO 2 H, -NHCH 2 C0 2 H, 
30 -C0 2 H, -CH 2 C0 2 H, -CH 2 CH 2 C0 2 H, -CH : CH 2 CH 2 C0 2 H, 
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-CH 2 CH(NH 2 )C0 2 H, -CH 2 CH[NHCOCH4> 2 ]C0 2 H, -CH 2 CH[NHCO(CH 2 ), 5 CH 3 
]C0 2 H, -CH 2 CH[NH-FMOC]C0 2 H, -CH 2 CH[NH-tBOC]C0 2 H, or a carboxylate. 
connected to the ring with a 0 to 3 carbon linker, 

5 -P0 3 H 2 , -CH 2 P0 3 H 2 , -CH 2 CH 2 P0 3 H 2 , -CH 2 CHNH 2 P0 3 H 2 . 

-CH 2 CH[NHCOCH^]P0 3 H 2 , -CH 2 CH[NHCO(CH 2 ) 15 CH 3 ]P0 3 H 2 , 
-CH 2 CH[NH-FMOC]P0 3 H 2 , -CH 2 CH[NH-tBOC]P0 3 H 2 , or a phosphate or 
phosphonate connected to the ring with a 0 to 3 carbon linker, 

10 -S0 3 H, -CH 2 S0 3 H, -CH 2 CH 2 S0 3 H, -CH 2 CHNH 2 S0 3 H, 

-CH 2 CH[NHCOCH<J> 2 ]S0 3 H, -CH 2 CH[NHCO(CH 2 ), 5 CH 3 ]S0 3 H, 
-CH 2 CH[NH-FMOC]S0 3 H, -CH 2 CH[NH-tBOC]S03H, or a sulfate or sulfite 
connected to the ring with a 0 to 3 carbon linker, 

15 or acts as the functional equivalent of CH,CH(NH 2 )C0 2 H of T3 in the molecular 

recognition domain when bound to a TR, wherein said R, can be optionally 
substituted with an amine, 

wherein R 2 is 

20 

H, halogen, CF 3 , OH, NH 2 , SH, CH 3 , -Et, 

or acts as the functional equivalent of H in the molecular recognition domain when 
bound to a TR, 

25 wherein R 3 is 

-H, halogen, -CF 3 , -OH, -NH 2 , -N 3 , -SH, -CH 3 , -Et, 

or acts as the functional equivalent of I in the molecular recognition domain when 
bound to a TR, 

30 

wherein R 5 is 
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92. 

-H, halogen, -CF 3 , -OH, -NH 2 , -N 3 , -SH, -CH 3 , -Et, or acts as the functional 
equivalent of I in the molecular recognition domain when bound to a TR, and R 3 
can be identical to R 5 , 

5 wherein R 6 is 



-H, halogen, -CF 3 , -OH, -NH 2 , -SH, -CH 3 , or acts as the functional equivalent of 
H in the molecular recognition domain when bound to a TR, and R 2 can be 
identical to R$, 

10 

wherein R 2 ' is 

-H, halogen, -CF 3 , -OH, -NH 2 , -N 3 , -SH, -CH 3 , -Et, or acts as the functional 
equivalent of H in the molecular recognition domain when bound to a TR, 

15 

wherein R 3 ' is any hydrophobic group, including 

halogen, -CF 3 , -SH, alkyl, aryl, 5- or 6-membered heterocyclic, cyano, or acts as 
the functional equivalent of I in the molecular recognition domain when bound to a 
20 TR, 

wherein R 4 ' is 

-H, halogen, -CF 3 , -OH, -NH 2 , NH 3 , -N(CH 3 ) 3 , carboxylate, phosphonate, 
25 phosphate or sulfate, -SH, -CH 3 , -Et, or akyl, aryl or 5- or 6-membered 

heterocyclic aromatic attached through urea or carbamate linkages to O or N or S 
at the R/ position, or acts as the functional equivalent of OH in the molecular 
recognition domain when bound to a TR, 

30 wherein R 5 ' is 
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-H, -OH, -NH 2 , -N(CH 3 ) 2 -SH -NH 3 , -N(CH 3 ) 3 . carboxylate, phosphonate, 
phosphate, sulfate, branched or straight chain alkyl having 1 to 9 carbons, 
substituted or unsubstituted aryl, wherein said substituted aryl is substituted with 
halogen or 1 to 5 carbon alkyl and wherein said aryl is optionally connected to the 
5 ring by a -CH r , aromatic heterocycle having 5 to 6 atoms, wherein said 

heterocycle may be substituted with one or more groups selected from -OH, -NH 2 , 
-SH, -NH 3 , -N(CH 3 ) 3 , carboxylate, phosphonate, phosphate, sulfate, heteroalkyl, 
heteroaryl, arylalkyl, heteroaryl alkyl, polyaromatic, polyheteroaromatic, wherein 
said R 5 ' may be substituted with polar or charged groups, 

10 

wherein R 6 ' is 

-H, halogen, -CF 3 , -OH, -NH 2 , -SH, -CH 3 , -Et, or acts as the functional 
equivalent of H in the molecular recognition domain when bound to a TR, 

15 

wherein X is 

O, S, S0 2 , NH, NR 7 , CH 2 , CHR 7 , CR 7 R 7 , wherein R 7 is alkyl, aryl or 5- or 
6-membered heterocyclic aromatic, 
20 and wherein said TR LBD ligand has an apparent Kd for binding TR LBD of 1 jtM or 
less. 

5. The crystal of claim 2, wherein said crystal of TR LBD ligand has the following 
unit cell dimensions in angstroms: a=117.00 ± 2%, b=80.00 ± 2%, c=63.00 ± 2%, 

25 at a 0 angle of 120.00 ± 2% degrees and a monoclinic space group C2. 

6. The crystal of claim 5, wherein said crystal coordinates said TR LBD ligand and 
said crystal is selected from the group consisting of crystals having one of the following 
sets of properties: 

30 
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1) unit cell dimensions in angstroms: a = 1 17. 16, b=80.52, c = 63.21, at a 0 
angle of 120.58 degrees, a monoclinic space group C2, 2.2 angstrom resolution 
and completeness value of 87.0%, 

5 2) unit cell dimensions in angstroms:. a = 117.19, b = 80.20, c = 63.23, at a 0 

angle of 120.60 degrees, a monoclinic space group C2, 2.0 angstrom resolution 
and completeness value of 82.4 %, and 

3) unit cell dimensions in angstroms: a = 117.18, b=80.12, c=63.13, at a j3 
10 angle of 120.69 degrees, a monoclinic space group C2, 2.2 angstrom resolution 

and completeness value of 93.7%. 

7. The crystal of claim 2, wherein said crystal further comprises a human protein. 

15 8. The crystal of claim 7, wherein said crystal has the following unit cell dimensions 
in angstroms: a = b=68.448 ± 2% and c= 130.559± 2% at an a angle of 90.00°, a 0 
angle of 90.00° and a y angle of 120.00° and having a trigonal space group p3(l)2L 

9. The crystal of claim 8, wherein said crystal further comprises a human protein. 
20 10. A computational method of designing a nuclear receptor synthetic ligand 
comprising: 

1) using a three dimensional model of a crystallized protein comprising a nuclear 
receptor LBD with a bound ligand to determine at least one interacting amino acid of the 
nuclear receptor LBD that interacts with at least one first chemical moiety of said bound 

25 ligand, and 

2) selecting at least one chemical modification of said first chemical moiety to 
produce a second chemical moiety with a structure that either decreases or increases an 
interaction between said interacting amino acid and said second chemical moiety 
compared to said interaction between said interacting amino acid and said first chemical 

30 moiety. 



11. The method of claim 10, wherein steps 1 and 2 are repeated. 
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12. The method of claim 10, further comprising generating said three dimensional 
model of said crystallized protein comprising a TR LBD with a bound TR ligand. 

13. The method of claim 12, wherein said three dimensional model is generated by 
5 comparing isomorphous ligand derivatives to produce improved phasing. 

14. The method of claim 13, wherein said isomorphous ligand derivatives are provided 
by substituting at least one of the R5, R3, 'R5 and 'R3 positions of a thyronine derivative 
with Br or I. 

10 

15. The method of claim 14, further comprising determining a change in interaction 
between said interacting amino acid and said ligand after chemical modification of said 
first chemical moiety. 

15 16. The method of claim 12, wherein said selecting uses said first chemical moiety 
that interacts with at least one of the said interacting amino acids listed in APPENDIX 2. 

17. The method of claim 16, wherein said chemical modification enhances hydrogen 
- bonding interaction, charge interaction, hydrophobic interaction, Van Der Waals 

20 interaction or dipole interaction between said second chemical moiety and said interacting 
amino acid compared to said first chemical moiety and said interacting amino acid. 

18. The method of claim 17, wherein said chemical modification reduces hydrogen 
bonding interaction, charge interaction, hydrophobic interaction, Van Der Waals 

25 interaction or dipole interaction between said second chemical moiety and said interacting 
amino acid compared to said first chemical moiety and said interacting amino acid. 



19. The method of claim 18, wherein said first chemical moiety is at least 4.5 
angstroms away from at least one distant amino acid and said distant amino acid is not 
30 any of the said interacting amino acids listed in claim 17. , 
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20. The method of claim 18, wherein said first chemical moiety is 6 to 12 angstroms 
away from a distant amino acid. 

21. The method of claim 19, wherein said chemical modification extends toward said 
5 distant amino acid and produces hydrogen bonding interaction, charge interaction, 

hydrophobic interaction, Van Der Waals interaction or dipole interaction between said 
second chemical moiety and said distant amino acid. 

22. The method of claim 19, wherein said chemical modification sterically hinders said 
10 distant amino acids from binding other amino acids but does not sterically hinder said 

interacting amino acids. 

23. The method of claim 17 or 19 further comprising determining a change in 
interaction between said interacting amino acid and said ligand after said chemical 

15 modification, wherein said determining and modifying are performed using a computer 
program to represent chemical structures of said interacting amino acid and ligand. 

24. The method of claim 22, wherein said chemical modification is at a R 5 ' position of 
a thyronine derivative. 

20 

25. The method of claim 24, wherein said chemical modification sterically hinders 
activation helix function. r 

26. The method of claim 25, wherein said chemical modification does not interfere 
25 with an interaction between said interacting amino acid and an atom from R 1 to 6 

positions of said ligand. 

27. The method of claim 26, wherein said activation helix is helix H12. 

30 28. The method of claim 27, wherein said chemical modification comprises a planar 
structure, equivalent in length to at least a 3 carbon alkyl that projects out from the plane 
of the prime ring of said thyronine derivative at least 30°. 
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29. The method of claim 10 further comprising generating said three dimensional 
model of a nuclear receptor other than TR using said nuclear receptor's LBD amino acid 
sequence and using said crystallized protein comprising a TR LBD with a bound TR LBD 
ligand. 

5 

30. The method of claim 29, wherein said chemical modification enhances hydrogen 
bonding interaction, charge interaction, hydrophobic interaction, Van Der Waals 
interaction or dipole interaction between said second chemical moiety and said interacting 
amino acid compared to said first chemical moiety and said interacting amino acid. 

10 

31. The method of claim 30, wherein said chemical modification reduces hydrogen 
bonding interaction, charge interaction, hydrophobic interaction, Van Der Waals 
interaction or dipole interaction between said second chemical moiety and said interacting 
amino acid compared to said first chemical moiety and said interacting amino acid. 

15 

32. The method of claim 31, wherein said chemical modification is only of said first 
chemical moiety that is at least 3 angstroms away from at least one distant amino acid and 
said distant amino acid is different from said interacting amino acid. 

20 33. The method of claim 33, wherein said first chemical moiety is 4 to 12 angstroms 
away from a distant amino acid. 

34. The method of claim 33, wherein said chemical modification extends toward said 
distant amino acid and produces hydrogen bonding interaction, charge interaction, 
25 hydrophobic interaction, Van Der Waals .interaction or dipole interaction between said 
second chemical moiety and said distant amino acid. 



30 



35. The method of claim 34, wherein said chemical modification sterically hinders said 
distant amino acids from binding other amino acids but does not sterically hinder said 
interacting amino acids. 
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36. The method of claim 35, wherein said first chemical moiety comprises C-H or 
C-OH. 

37. The method of claim 36, wherein said second chemical moiety comprises a longest 
entity equivalent to at least 4 double bonded carbons in length and a volume of at least 4 

5 double bonded carbons. 

38. A computational method of designing a nuclear receptor antagonist from a nuclear 
receptor agonist comprising: 

1) determining a structure of a molecular recognition domain of said agonist using 
10 a three dimensional model of a crystallized protein comprising a nuclear receptor LBD, 

and 

2) selecting at least one chemical modification of said agonist that provides a 
ligand structure that extends beyond a binding site for said agonist and in the direction of 
at least one protein domain important in nuclear receptor function. 

15 

39. The method of claim 38, wherein said protein domain is either: 

a) a transcription activation domain of said LBD, 

b) a repressor binding domain of said LBD, 

c) a DNA binding domain of said nuclear receptor, 

20 d) a heat shock protein binding domain of said nuclear receptor, 

e) a dimerization domain of said LBD, or 

f) a hinge region to said DNA binding domain. 

40. The method of claim 38, wherein said crystallized protein comprises said nuclear 
25 receptor LBD bound to a nuclear receptor ligand. 

41. The method of claim 40 further comprising generating said three dimensional 
model of said crystallized protein comprising a TR LBD with a bound TR ligand. 



30 42. 



The method of claim 40 wherein said three dimensional model is for a TR LBD. 
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43. The method of claim 41 further comprising generating said three dimensional 
model of a nuclear receptor other than TR using said nuclear receptor's LBD amino acid 
sequence. 

5 44. The method of claim 39 wherein said three dimensional model is for an ER LBD. 

45. The method of claim 42 wherein said LBD is from a receptor selected from a 
group consisting of glucocorticoid receptor, estrogen receptor, retinoid receptor and 
vitamin D receptor. 

10 

46. The method of claim 39, wherein said chemical modification minimizes a loss of 
hydrogen bonding interaction, charge interaction, hydrophobic interaction. Van Der 
Waals interaction or dipole interaction between said molecular recognition domain and 
said antagonist compared agonist hydrogen bonding interaction, charge interaction, 

15 hydrophobic interaction or dipole interaction with said molecular recognition domain. 

47. The method of claim 39, wherein said chemical modification reduces hydrogen 
bonding interaction, charge interaction, hydrophobic interaction or dipole interaction 
between said molecular recognition domain and said antagonist compared with the agonist 

20 hydrogen bonding interaction, charge interaction, hydrophobic interaction. Van Der 
Waals interaction or dipole interaction with said molecular recognition domain but said 
chemical modification still permits said antagonist to bind to said nuclear receptor with an 
affinity of lOOnM Kd or less. 

25 48. The method of claim 39 or 47, wherein said molecular recognition domain is at 
least 3 angstroms away from at least one distant amino acid and said distant amino acid 
does not contribute significantly to binding said molecular recognition domain. 



49. The method of claim 48, wherein said first molecular recognition domain is 4 to 
30 12 angstroms away from a distant amino acid. 
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50. The method of claim 49, wherein said chemical modification extends toward said 
distant amino acid and produces a hydrogen bonding interaction, charge interaction, 
hydrophobic interaction, Van Der Waals interaction or dipole interaction between said 
distant amino acid and said molecular recognition domain. 

5 

51. The method of claim 50, wherein said chemical modification sterically hinders said 
distant amino acids from binding other amino acids but does not sterically hinder amino 
acids that bind said molecular recognition domain. 

10 52. The method of claim 51, wherein said chemical modification of said ligand is at a 
C-H or C-OH. 

53. The method of claim 52, wherein said second chemical modification consists of a 
longest entity equivalent to at. least 4 double bonded carbons in length and has volume of 

15 at least 4 double bonded carbons in length. 

54. The method of claim 44, wherein said three dimensional model is represented in 
two dimensions. 

20 55, A computational method of designing a nuclear receptor super agonist or 
antagonist comprising: 

1) determining at least one interacting amino acid of a nuclear receptor LBD that 
interacts with at least one first chemical moiety of said ligand using a three dimensional 
model of a crystallized protein comprising a nuclear receptor LBD with a bound ligand, 

25 and 

2) selecting at least one chemical modification of said first chemical moiety to 
produce a second chemical moiety with a structure to reduce or enhance an interaction 
between said interacting amino acid and said second chemical moiety compared to said 
interaction between said interacting amino acid and said first chemical moiety. 

30 



56. The method of claim 55, wherein said chemical modification enhances hydrogen 
bonding interaction, electrostatic interaction, charge interaction, hydrophobic interaction. 
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said interacting amino acid compared to said first chemical moiety and said interacting 
amino acid. 

5 57. The method of claim 55, wherein said chemical modification changes a 

carboxylate moiety of said first chemical moiety to a phosphonate or phosphate to make 
said second chemical moiety. 

58. The method of claim 55, wherein said nuclear receptor is TR and said chemical 
10 modification enhances said interaction between said second chemical moiety and at least 
one of the following arginines: Arg 262, Arg 266 or Arg 228 of the rat a-TR or an 
arginine of human a-TR or 0-TR that corresponds in its three dimensional position in said 
three dimensional model to either said arginines: Arg 262, Arg 266 or Arg 228. 

15 59. The method of claim 58, wherein said chemical modification changes a side group 
to fill space normally occupied by water when an agonist or naturally occurring ligand is 
bound to said LBD. 

60. The method of claim 59, wherein said chemical modification snugly fits into a 
20 pocket or surface of a binding cavity of said LBD and complements to a charge or 

hydrophobicity or both of said pocket or surface. 

61. The method of claim 60, wherein said chemical modification changes said first 
chemical moiety either a first carbon bonded either to 1) a first phenyl at Rl of T3 or 2) 

25 a first phenyl at Rl of a T3 agonist with a T3 carboxylate group to a second chemical 
moiety a two carbon space group. 

62. The method of claim 55, wherein said chemical modification changes said first 
chemical moiety either 1) a carbon linking a first phenyl at R4 of T3 to a second phenyl 

30 ring of T3 or 2) a carbon linking a first phenyl at R4 of a T3 agonist with a second 
phenyl ring of said T3 agonist to a second chemical moiety mono- or geminal-substituted 
carbon group. 
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63. The method of claim 55, wherein said chemical modification changes said first 
chemical moiety either 1) a carbon linking a first phenyl at R4 of T3 to a second phenyl 
ring of T3 or 2) a carbon linking a first phenyl at R4 of a T3 agonist with a second 
phenyl ring of said T3 agonist. 
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64. A compound of the formula: 
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103- 



10 



wherein R, is 
15 -0-CH 2 CO 2 H, -NHCH 2 C0 2 H, 

-C0 2 H, -CH 2 C0 2 H, -CH 2 CH 2 C0 2 H, -CH 2 CH 2 CH 2 C0 2 H, 
-CH 2 CH(NH 2 )C0 2 H, -CH 2 CH[NHCOCH<f> 2 ]CO,H, -CH : CH[NHCO(CH 2 ) 15 CH 
]C0 2 H, -CH 2 CH[NH-FMOC]C0 2 H, -CH 2 CH[NH-tBOC]C0 2 H, or a carboxylate 
connected to the ring with a 0 to 3 carbon linker, 



3 



20 



25 



-P0 3 H 2 , -CH 2 P0 3 H 2) -CH 2 CH 2 P0 3 H 2 , -CH 2 CHNH 2 P0 3 H 2 , 
-CH 2 CH[NHCOCH4> 2 ]P0 3 H 2 , -CH 2 CH[NHCO(CH 2 ) 15 CH 3 ]P0 3 H 2 , 
-CH 2 CH[NH-FMOC]P0 3 H 2 , -CH 2 CH[NH-tBOC]P0 3 H 2 , or a phosphate or 
phosphonate connected to the ring with a 0 to 3 carbon linker, 

-S0 3 H, -CH 2 S0 3 H, -CH 2 CH 2 S0 3 H, -CH 2 CHNH 2 S0 3 H, 
-CH 2 CH[NHCOCH«fc]S0 3 H, -CH 2 CH[NHCO(CH 2 ) l5 CH 3 ]S0 3 H, 
-CH 2 CH[NH-FMOC]S0 3 H, -CH 2 CH[NH-tBOC]S03H, or a sulfate or sulfite 
connected to the ring with a 0 to 3 carbon linker, 



30 
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or acts as the functional equivalent of CH 2 CH(NH 2 )C0 2 H of T3 in the molecular 
recognition domain when bound to a TR, wherein said R, can be optionally 
substituted with an amine, 

5 wherein R 2 is 

H, halogen, CF 3 , OH, NH 2 , SH, CH 3 , -Et, 

or acts as the functional equivalent of H in the molecular recognition domain when 
bound to a TR, 

10 

wherein R 3 is 

-H, halogen, -CF 3 , -OH, -NH 2 , -N 3 , -SH, -CH 3 , -Et, 

or acts as the functional equivalent of I in the molecular recognition domain when 
15 bound to a TR, 

wherein R 5 is 

-H, halogen, -CF 3 , -OH, -NH 2 , -N 3 , -SH, -CH 3 , -Et, or acts as the functional 
20 equivalent of I in the molecular recognition domain when bound to a TR, and R 3 

can be identical to R 5 , 

wherein R 6 is 

25 -H, halogen, -CF 3 , -OH, -NH 2 , -SH, -CH 3 , or acts as the functional equivalent of 

H in the molecular recognition domain when bound to a TR, and R 2 can be 
identical to R$, 

wherein R 2 ' is 

30 

-H, halogen, -CF 3 , -OH, -NH 2 , -N 3 , -SH, -CH 3 , -Et, or acts as the functional 
equivalent of H in the molecular recognition domain when bound to a TR, 
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wherein R 3 ' is any hydrophobic group, including 

halogen, -CF 3 , -SH, alkyl, aryl, 5- or 6-membered heterocyclie, cyano, or acts as 
the functional equivalent of I in the molecular recognition domain when bound to a 
5 TR, 

wherein R 4 ' is 

-H, halogen. -CF 3 , -OH, -NH 2 , NH 3 , -N(CH 3 ) 3 , carboxylate, phosphonate, 
10 ■ phosphate or sulfate, -SH, -CH 3 , -Et, or akyl, aryl or 5- or 6-membered 

heterocyclic aromatic attached through urea or carbamate linkages to 0 or N or S 
at the R/ position, or acts as the functional equivalent of OH in the molecular 
recognition domain when bound to a TR, 

15 wherein R 5 ' is 

-H, -OH, -NH 2 , -N(CH 3 ) 2 -SH -NH 3 , -N(CH 3 ) 3 , carboxylate, phosphonate, 
phosphate, sulfate, branched or straight chain alkyl having 1 to 9 carbons, 
substituted or unsubstituted aryl, wherein said substituted aryl is substituted with 

20 halogen or 1 to 5 carbon alkyl and wherein said aryl is optionally connected to the 

ring by a -CH 2 -, aromatic heterocycle having 5 to 6 atoms, wherein said 
heterocycle may be substituted with one or more groups selected from -OH, -NH,, 
-SH, -NH 3 , -N(CH 3 ) 3 , carboxylate, phosphonate, phosphate or sulfate, heteroalkyl, 
arylalkyl, heteroaryl alkyl, polyaromatic, or polyheteroaromatic, wherein said R s ' 

25 may be substituted with polar or charged groups, 

wherein R 6 ' is 

-H, halogen, -CFj, -OH, -NH 2 , -SH, -CH 3 , -Et, or acts as the functional 
30 equivalent of H in the molecular recognition domain when bound to a TR, 



wherein X is 
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O, S, S0 2 , NH, NR 7l CH 2 , CHR ? , CR 7 R 7l wherein R 7 is alkyl, aryl or 5- or 
6-membered heterocyclic aromatic, 

and wherein said TR LBD ligand has an apparent Kd for binding TR LBD of 1 or 
5 less. 

65. The compound of claim 64, wherein 

R, is carboxylate, phosphonate, phosphate or sulfite and is connected to the ring with a 
0 to 3 carbon linker, 
10 R 2 is H, 

R 3 is -I, -Br, or -CH 3 , 
R 5 is -I, -Br, or -CH 3 , 
R6is -H, 
R 2 ' ^ -H, 

15 R 3 ' is -I, -Br, -CH 3 , -iPr, -phenyl, benzyl, or 5- or 6-membered ring heterocycles, 
R/ is -OH, -NH 2f and -SH, 

R 5 ' is -H, -OH, -NH 2 , -N(CH 3 ) 2 -SH -NH 3 , -N(CH 3 ) 3 , carboxylate, phosphonate, 
phosphate, sulfate, branched or straight chain alkyl having 1 to 9 carbons, 
substituted or unsubstituted aryl, wherein said substituted aryl is substituted with 

20 halogen or 1 to 5 carbon alkyl and wherein said aryl is optionally connected to the 

ring by a -CH 2 -, aromatic heterocycle having 5 to 6 atoms, wherein said 
heterocycle may be substituted with one or more groups selected from -OH, -NH 2 , 
-SH, -NH 3 , -N(CH 3 ) 3 , carboxylate, phosphonate, phosphate or sulfate, heteroaikyl, 
arylalkyl, heteroaryl alkyl, polyaromatic, or polyheteroaromatic, wherein said R 5 ' 

25 may be substituted with polar or charged groups, and 

R*' is H. 

66. The compound of claim 65, wherein said compound is made by the method of 
claim 13. 

30 67. The compound of claim 65, wherein said compound is made by the method of 
claim 39. 

68. The compound of claim 65, wherein said compound is a TR antagonist. 
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69. The compound of claim 65, wherein said compound is a TR agonist. 

70. The compound of claim 65, wherein said compound is a TR a selective ligand. 

71. The compound of claim 65, wherein said compound is a TR j3 selective ligand. 

72. A pharmaceutical composition having selective thyromimetic activity comprising a 
5 compound as claimed in claim 65 and a pharmaceutical^ effective carrier. 

73. The composition of claim 72, wherein R5' is an alkyl having from 1 to 9 carbons 
and being straight chain or branched. 

74. A method of lowering the ratio of LDL-cholesterol to HDL-cholesterol levels 
comprising administering an effective amount of the compound of claim 65. 

10 75. A method of lowering plasma lipid levels which comprises administering to an 
animal in need thereof an effective amount of the compound of claim 65. 
76. A method of treating thyroid hormone deficiency in patients with compromised 
cardiac function, comprising administering an effective amount of the compound of claim 
65. 
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RECEPTORS 

Figure shows steps in the actions of the nuclear hormone receptors. Ligands 
for these receptors circulate bound to varying extents to plasma proteins. 
Following dissociation from these proteins, these ligands enter ceils and bind 
specifically to their receptors. Steroids and vitamin D probably enter cells 
through passive di/fusion; thyroid he cm one and rednoic acid entry might 
involve specific transport processes (5). 
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